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^ ■ ABSTRACT 

' The bright Seyfert 1 galaxy MCG-6-30-15 has provided some of the best evidence to date 

for the existence of supermassive black holes in active galactic nuclei. Observations with 
ASCA revealed an X-ray iron line profile shaped by strong Doppler and gravitational effects. 
In this paper the shape of the iron line, its variability characteristics and the robustness of this 
spectral interpretation are examined using the long XMM-Newton observation taken in 2001. 
A variety of spectral models, both including and excluding the effects of strong gravity, are 
^-j. [ compared to the data in a uniform fashion. The results strongly favour models in which the 

■ spectrum is shaped by emission from a relativistic accretion disc. It is far more difficult to 
^-H ' explain the 3 — lOkeV spectrum using models dominated by absorption (either by warm or 
CO . partially covering cold matter), emission line blends, curved continua or additional contin- 

■ uum components. These provide a substantially worse fit to the data and fail to explain other 

observations (such as the simultaneous BeppoSAX spectrum). This reaffirms the veracity of 

Q^. the relativistic 'disc line' interpretation. The short term variability in the shape of the energy 

I ' spectrum is investigated and explained in terms of a two-component emission model. Using 

O \ a combination of spectral variability analyses the spectrum is successfully decomposed into a 

. variable power-law component (PLC) and a reflection dominated component (RDC). The for- 

^ ' mer is highly variable while the latter is approximately constant throughout the observation, 

. . , leading to the well-known spectral variability patterns. Consideration of the long term X-ray 

■ monitoring of MCG-6-30-15 by RXTE demonstrates that the long XMM-Newton observation 
L>( I sampled the 'typical' state of the source. These results and those of other analyses of the large 

. XMM-Newton dataset are summarised and their implications for understanding this enigmatic 

■ source are discussed. 
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1 INTRODUCTION 

X-ray spectroscopy holds great promise for probing physics close 
to the accreting black holes thought to power Active Galactic Nu- 
clei (AGN) and Galactic Black Hole Candidates (GBHCs). The 
high X-ray luminosities observed from these systems are thought 
to be generated close to the black hole (<; 100 rg; where rg = 
GAl/(? is the gravitational radius) where the effects of strong grav- 
ity, such as gravitational redshift and light bending, will imprint 
characteristic signatures on the emerging X-ray spectrum. 

As it escapes towards the observer, the primary X-ray emis- 
sion will interact with material in the innermost regions, modify- 
ing the observed X-ray spectrum. Reprocessing of the X-rays by 
an accretion disc will lead to a particular set of spectral features. 
Absorption, fluorescence, recombination, and Compton scattering 
in the surface layers of the X-ray illuminated disc produce a 're- 
flection' spectrum. Under a wide range of physical conditions the 



most prominent observables will be an emission line at « 6.4 keV 
(the iron Kq fluorescence line) and a broad 'hump' in the spectrum 
which peaks in the range 20 — 40 keV (e.g. Guilbert & Rees 1988; 
Lightman & White 1988; George & Fabian 1991; Matt, Perola & 
Piro 1991). Observations of Seyfert 1 galaxies with Ginga (Pounds 
et al. 1990; Nandra & Pounds 1994) demonstrated the simultane- 
ous presence of a ~ 6.4 keV emission line and an up-turn in the 
spectrum above ^ 10 keV, indicative of the reflection hump. The 
presence of the reflection hump has since been confirmed by CGRO 
(e.g. Zdziarski et al. 1995) and BeppoSAX (e.g. Perola et al. 2002). 

The bright Seyfert 1 galaxy MCG-6-30-15 {z = 0.007749) 
has received particular attention since ASCA observations revealed 
the presence of a broad, asymmetric emission feature at energies 
between 4 — 7 keV identified with a highly broadened Fe Kq emis- 
sion line (Tanaka et al. 1995). Similar features were subsequently 
observed in other Seyfert 1 galaxies (Mushotzky et al. 1995; Nan- 
dra et al. 1997). The profile of the line can be explained in terms 
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of fluorescent emission from the surface of an accretion disc ex- 
tending down to < 6 Tg about a supermassive black hole (SMBC): 
the relativistic 'disc line' model (Fabian et al. 1989; Stella 1990; 
Laor 1991). For comparison, the emission lines seen in other wave- 
bands (e.g. broad optical lines) are thought to originate at distances 
> 10^ r-g. The broad Fe Ka line therefore potentially offers a pow- 
erful diagnostic of the physical conditions in the immediate envi- 
ronment of the black hole (see Fabian et al. 2000; Nandra 2001 and 
Reynolds & Nowak 2003 for recent reviews). 

Repeated observations of MCG-6-30-15 withA5CA (Iwasawa 
et al. 1996, 1999; Shih et al. 2002), BeppoSAX (Guainazzi et al. 
1999), RXTE (Lee et al. 1999, 2000; Vaughan & Edelson 2001), 
Chandra (Lee et al. 2002) and XMM-Newton (Wilms et al. 2001, 
hereafter WOl, and Fabian et al. 2002, hereafter paper I) have 
confirmed the presence of the broad, asymmetric emission fea- 
ture. However, while this is clearly not an artifact of any partic- 
ular instrument or observation, there are uncertainties associated 
with the detailed modelling of the spectrum. The most challenging 
of these is that MCG-6-30-15, along with many other Seyfert Is, 
shows complex absorption in its X-ray spectrum. This complicates 
the process of identifying the correct underlying continuum and 
thereby measuring the superposed line emission. Indeed, the pres- 
ence of relativistically broadened emission lines in the X-ray spec- 
tra of other Seyfert 1 galaxies has recently become something of a 
cause celebre (see e.g. Lubinski & Zdziarski 2001; Inoue & Mat- 
sumoto 2001; Branduardi-Raymont et al. 2001; Pounds & Reeves 
2002; Page, Davis & Salvi 2003; Pounds et al. 2003a,b). 

The long (~ 320 ks) XMM-Newton observation of MCG~6- 
30-15 is of great importance to studies of the broad iron line, since 
MCG-6-30-15 offers the best established example of a relativistic 
line profile. In this paper the X-ray spectrum and spectral variabil- 
ity of MCG-6-30-15 are examined using the long XMM-Newton 
observation. The interpretation of the spectrum in terms of rela- 
tivistic Fe Ka emission is investigated by fitting a range of models 
including both simple phenomenological models and models ac- 
counting for the detailed physics expected in photoionised accre- 
tion discs. These are compared to models that do not include rela- 
tivistic effects in order to test the strength of the argument for strong 
gravity in MCG-6-30-15. A variety of techniques are employed to 
characterise the spectral variability of the source and thereby fur- 
ther constrain the range of viable spectral models. The two XMM- 
Newton observations are considered in the context of the long term 
behaviour of MCG~6-30-I5 using RXTE monitoring observations 
spanning more than six years. 

The rest of this paper is organised as follows. Section|2|gives 
details of the XMM-Newton and RXTE observations and their data 
reduction. Section|3|briefly reviews the RXTE monitoring of MCG- 
6-30-15 and shows how the two XMM-Newton observations fit in 
with the long timescale properties of the source. Section|4]discusses 
the spectral fitting results using models that include relativistic ef- 
fects (Section f4.6> and exclude relativistic effects (Section [4.7t . In 
section |5| the spectral variability of the source is examined using 
flux-flux plots (section [5. 1> . flux-resolved spectral analysis (sec- 
tion rms spectra (section f5.3> . time-resolved spectral anal- 
ysis (section [5 .4t and Principal Component Analysis (PCA; sec- 
tion l5.5> . Section|6|discusses differences between the long XMM- 
Newton observation of 2001 and the previous observation taken in 
2000. Section|7|summarises the new results from the present work 
as well as the other analyses of these XMM-Newton data. Section|8| 
presents a discussion of the main results and the main conclusions 
are given in section|9| Some details of the EPIC spectral calibration 
are investigated in Appendix IaI 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 XMM-Newton observations 

XMM-Newton (Jansen et al. 2001) has observed MCG-6-30-15 
twice. The first observation was performed during revolution 108 
(2000 July 11-12) and the second observation was performed dur- 
ing revolutions 301-303 (2001 July 31 - August 5). In both cases 
the source was observed on-axis and all instruments were operat- 
ing nominally. The first results from the 2000 XMM-Newton obser- 
vation were presented in WOl and the first results from the 2001 
XMM-Newton observation were presented in paper I. 

During the 2000 observation the pn camera was operated in 
small window mode, the MOS2 camera in full frame mode and the 
MOSI camera in timing mode. During the 2001 observation both 
the EPIC MOS cameras and the EPIC pn camera were operated in 
small window mode. For the present analysis only the EPIC data 
taken in small window mode were used (i.e. only the pn from the 
2000 observation, pn and both MOS from the 2001 observation). 

The pn small window mode uses a 63 x 64 pixel window 
(~ 4'. 3 X 4'. 4 on the sky) with the source positioned less than an 
arcminute from the CCD boundary. The MOS small window mode 
uses a window of 100 x 100 pixels (~ 1'.8 x 1'.8) with the source 
image approximately centred within the window. The small win- 
dow modes allow CCD frames to be read out every 5.7 ms for the 
pn and 0.3 s for the MOS. These short frame times help lessen the 
effects of photon pile-up which can distort the spectrum of bright 
sources (Ballet 1999; Lumb 2000). All three EPIC instruments used 
the medium filter. 



2.2 XMM-Newton data reduction 

The extraction of science products from the Observation Data Files 
(ODFs) followed standard procedures using the XMM-Newton Sci- 
ence Analysis System (SAS v5 . 4 . 1). The EPIC data were pro- 
cessed using the standard SAS processing chains to produce cal- 
ibrated event lists. These removed events from the position of 
known defective pixels, corrected for Charge Transfer Inefficiency 
(CTI) and applied a gain calibration to produce a photon energy for 
each event. 

Source data were extracted from circular regions of radius 
40" from the processed images and background events were ex- 
tracted from regions in the small window least effected by source 
photons. These showed the background to be relatively low and 
stable throughout the observations, with the exception of the final 
few ks of each revolution where the background rate increased (as 
the spacecraft approached the radiation belts at perigee). Data from 
these periods were ignored. The total amount of 'good' exposure 
time selected was 64 ks from the pn during the 2000 observation, 
and for the 2001 observation 315 ks and 228 ks from the MOS 
and pn, respectively. (The lower pn exposure is due to the lower 
'live time' of the pn camera in small-window mode, ~ 71 per cent; 
Striideret al. 2001). 

The ratios of event patterns as a function of energy showed 
there is negligible pile-up in the pn data but the MOS data suffer 
slightly from pile-up. This may lead to a slight distortion of the 
MOS spectra. In order to minimise this effect only photons cor- 
responding to single pixel events (EPIC patterns 0) were extracted 
from the MOS cameras. The spectrum of single pixel events should 
be least affected by pile-up (Molendi & Sembay 2003). For the pn, 
single pixel (pattern 0) and double pixel (patterns 1 — 4) events 
were extracted separately for the spectral analysis and fitted simul- 
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Figure 1. Top panel: the 2 — lOkcV RXTE light curve around the time of the two XMM-Newton observations (during PCA gain epoch 5). Bottom panel: 
Hardness ratio (2 — 4/4 — lOkeV). The solid line marks the weighted mean hardness ratio throughout the observing period. 



taneously (the spectrum of single pixel events provides better en- 
ergy resolution). Response matrices were generated using RMFGEN 
vl . 48 . 5 and ancillary responses were generated with ARFGEN 
V . 1 . 5 4 . 7 . For the 2001 observation source spectra were grouped 
such that each energy bin contains at least 2, 000 counts for the pn 
and 400 counts for the MOS^. For the shorter 2000 observation the 
pn spectrum was grouped to contain at least 500 counts per energy 
bin. Background subtracted light curves were extracted from the pn 
camera in 1000 s bins using both single and double pixel events. 

2.3 RXTE observations 

MCG-6-30-15 has been regularly monitored by RXTE since 1996. 
The RXTE Proportional Counter Array (PCA) consists of five col- 
limated Proportional Counter Units (PCUs), sensitive to X-rays in 
a nominal 2 — 60 keV bandpass and with a total collecting area of 
~ 6250 cm^ . However, due to operational constraints there are usu- 
ally fewer than five PCUs operated during any given observation. 
Only data from PCU2 were analysed here as this is the one unit that 
has been operated consistently throughout the lifetime of the mis- 
sion. In the following analysis, data from the top (most sensitive) 
layer of the PCU array were extracted using SAEXTRCT v4 . 2d. 
Poor quality data were excluded using standard selection criteria^. 
The PCA background was estimated using the latest version of the 

^ This grouping was used instead of the more traditional A'^ = 20 as it 
results in fewer spectral bins (and therefore reduces the computation time 
for each fit iteration) but does not compromise the intrinsic spectral reso- 
lution. The data contain sufficient counts that even with this heavy binning 
the instrument resolution is oversampled. 

^ The acceptance criteria were as follows: TIME_SINCE_SAA> 20 min; 
Earth elevation angle ELV!^ 10°; offset from optical position of MCG- 



'combined' model (Jahoda et al. 2000). Light curves were initially 
extracted from the STANDARD-2 data with 16 s time resolution. 
The data were rebinned into quasi-continuous observations lasting 
typically 1 — 2 ks. 



3 X-RAY HISTORY OF MCG-6-30-15 

FigureQshows the RXTE light curve and hardness ratio around the 
times of the two XMM-Newton observations. From this it is clear 
that the source displays strong X-ray variability, but most is con- 
fined to short term 'flickering' rather than longer term trends in the 
data. A counter example can be seen around MJD 52,000 when 
there appears to have been a prolonged period of lower flux and 
variability, during which the spectrum hardened slightly. However, 
both XMM-Newton observations seem to have caught MCG-6-30- 
15 in its 'typical' state; there are no outstanding secular changes in 
either flux or hardness ratio to suggest anything other than normal 
behaviour in the source during the two XMM-Newton observations. 
The same conclusion is reached when the RXTE monitoring ex- 
tending back to 1996 are considered. Papadakis et al. (2002) and 
Uttley, M'^Hardy & Papadakis (2002) discuss the RXTE monitoring 
observations in more detail. 



4 SPECTRAL FITTING ANALYSIS 

This section describes the results of the spectral fitting analysis. 
The raw spectra from the 2001 observation are plotted in Fig.|2| A 

6-30-15 OFFSET < 0.02°; and ELECTRON2 ^ 0.1. This last criterion re- 
moves data with high anti-coincidence rate in the propane layer of the PCU. 
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Figure 2. Full-band EPIC spectra. The four sets of data are the pn spectra for single and double events (shown separately) and the MOSl and M0S2 spectra. 
For clarity the MOS2 spectrum has been shifted downwards by 20 per cent. Marked are the positions of the most significant instrumental features (edges due 
to neutral O, Si and Au) and the positions of the two features expected to be strong in the source spectrum (the O VII edge and Fe Ka emission line). 



variety of models, both including and excluding relativistic effects, 
were fitted to the data. The goal of the analysis was to find the 
models that can accurately match the data, and test whether models 
including relativistic effects were preferred. This analysis expands 
upon the analyses of the same observation presented in paper I, 
Fabian & Vaughan (2003) and Ballantyne et al. (2003) in that the 
latest software and calibration were used (as of 2003 May), a wider 
variety of spectral models was tested, and more information from 
other observations was used to better constrain the range of the 
model fits. 



4.1 Limits of spectral modelling 

The high signal-to-noise of these spectra mean that systematic er- 
rors in the detector response model can dominate over the statistical 
errors (due to photon noise). A brief analysis of some calibration 
data is discussed in Appendix]^ This demonstrated that signifi- 
cant instrumental features should be expected close to the K-edges 
of O and Si and the M-edge of Au. Elsewhere the spectrum shows 
no sharp instrumental residuals larger than ~ 5 per cent. Therefore 
in the spectral fitting outlined below, the fitting is primarily carried 
out over the 3 — 10 keV spectral range that excludes these features. 
Appendix IaI also demonstrates that there is a strong discrepancy 
between the spectral slopes derived from each of the EPIC cameras 
(see Molendi & Sembay 2003). Therefore the spectral slope of the 
power-law continuum was always allowed to vary independently 
between the three EPIC cameras. 

With these precautions in place various spectral models were 
compared to the data as described below. The models were fitted 
to the four spectra (pn:singles, pn:doubles, MOSl and MOS2) si- 
multaneously using XSPEC vll.2 (Arnaud 1996). The quoted 



errors on the derived model parameters correspond to a 90 per cent 
confidence level for one interesting parameter (i.e. a A^^ = 2.71 
criterion), unless otherwise stated, and fit parameters are quoted for 
the rest frame of the source. 



4.2 A first look at the spectrum 

The count spectra shown in Fig. |2| give only a poor indication of 
the true spectrum of the source since they are strongly distorted 
by the response of the detector. In order to gain some insight into 
the broad-band shape of the source spectrum, without recourse to 
model fitting, the EPIC data for MCG~6-30-15 were compared to 
the spectrum of 3C 273 (see Appendix]^. The raw EPIC pn:s data 
(shown in Fig.|2j were divided by the raw pn:s spectrum for 3C 273. 
This gave the ratio of the spectra for the two sources. 3C 273 was 
chosen as it is a bright source and has a relatively simple spectrum 
in the EPIC band, i.e. a hard power-law plus smooth soft excess 
modified by Galactic absorption (Page et al. 2003b), and in par- 
ticular does not contain any strong, sharp spectral features such as 
lines or edges. Thus the ratio of MCG-6-30-15 to 3C 273 spectra 
will factor out the broad instrumental response to give a better im- 
pression of the true shape of the MCG~6-30-15 spectrum. The top 
panel of Fig. |3| shows this ratio. The advantage of this plot is that 
it gives an impression of the MCG-6-30-15 spectrum and is com- 
pletely independent of spectral fitting. The disadvantages are that it 
is is in arbitrary units and only gives a very crude representation of 
the true flux spectrum of MCG-6-30-15 because instead of being 
distorted by the instrumental response, it is now distorted by the 
spectrum of 3C 273. 

A more accurate, but model-dependent method for obtain- 
ing a 'fluxed' spectrum of MCG-6-30-15 is to normalise the ra- 
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Figure 3. EPIC pn:s spectrum of MCG-6-30-15 after normalising by the 
3C 273 spectrum. The top panel shows the ratio of counts in the MCG-6-30- 
15 and 3C 273 spectra. The bottom panel shows the ratio normalised by the 
spectral model for 3C 273. This reveals the 'fluxed' spectrum of MCG-6- 
30-15 (in units of keV^ cm^ keV~^). The 'fluxed' spectrum has also 
been corrected for Galactic absorption. (The sharp feature at fs 0.54 keV 
is an ailifact caused by the neutral O K-edge in the absorption model being 
sharper than the instrumental response.) 



tio by a spectral model for 3C 273. This technique is very similar 
to the method routinely used in optical spectroscopy of normalis- 
ing the target spectrum using a standard star spectrum to obtain a 
wavelength-dependent flux calibration. The MCG~6-30-15/3C 273 
ratio was multiplied by a spectral model for 3C 273 (defined in 
flux units). The spectral model comprised a hard power-law plus 
two blackbodies to model the soft excess, modified by Galactic ab- 
sorption (see AppendixlAland also Page et al. 2003b). This process 
transformed the raw MCG~6-30-15 count spectrum into a 'fluxed' 
spectrum, without directly fitting the MCG~6-30-15 data. (The ac- 
curacy of this procedure does however depend on the accuracy of 
the 3C 273 spectral modelling.) The fluxed spectrum for MCG~6- 
30-15 is shown in the bottom panel of Fig.|3| 

Note that this procedure is not a formal deconvolution of the 
spectrum (see Blissett & Cruise 1979; Kahn & Blissett 1980) and 
so does not correctly take into account the finite detector spectral 
resolution. Nevertheless, it does provide a simple yet revealing im- 
pression of the source spectrum. Specifically, the plot clearly re- 
veals the strongest spectral features in MCG-6-30-15, namely the 
strong warm absorption concentrated in the ~ 0.7 — 2 keV region 
and the iron line peaking at « 6.4 keV. 



4.3 Ingredients for a realistic spectral model 

4.3.1 Soft X-ray absorption 

In the spectral fitting analysis described below the absorption due to 
neutral gas along the line-of-sight through the Galaxy is accounted 
for assuming a column density of A^h = 4.06 x 10^" cm~^ 
(derived from 21cm measurements; Elvis, Wilkes & Lockman 
1989). Accurately modelling the intrinsic absorption provided a 
more challenging problem. MCG-6-30-15 is known to possess a 
complex soft X-ray spectrum, affected by absorption in partially 
ionised ('warm') material (Nandra & Pounds 1992; Reynolds et al. 
1995; Otani et al. 1996) and possibly also dust (Reynolds et al. 
1997; Lee et al. 2001; Turner et al. 2003a; Ballantyne, Weingart- 



ner & Murray 2003). However, the detailed structure of the soft 
X-ray spectrum, as seen at improved spectral resolution with the 
grating instruments on board Chandra and XMM-Newton, remains 
controversial (Branduardi-Raymont 2001; Lee et al. 2001; Sako et 
al. 2003; Turner et al. 2003a). In the present analysis the complicat- 
ing effects of the warm absorber were mitigated by concentrating 
on the spectrum above 3 keV. However, this will not entirely elim- 
inate the effects of absorption. 

In the 3 — 10 keV range absorption hy Z ^ 14 ions will be 
limited since their photoionisation threshold energies are less than 
3keV for all charge states (the ionisation energy for H-like Si^^'^ 
is 2.67 keV; Verner & Yakovlev 1995). Similarly, the L-shell of 
ionised Fe can significantly absorb the spectrum around 1 keV due 
to complexes of absorption lines and edges (e.g. Nicastro, Fiore & 
Matt 1999; Behar, Sako & Kahn 2001). These transitions all oc- 
cur below 2.05 keV (the ionisation energy of the L-edge in Li-like 
Fe^^''). Thus above 3keV the only opacity from \ow-Z ions and 
Fe L will be due to the roll-over of their superposed photoelectric 
absorption edges. The dominant effect of the warm absorber on the 
spectrum above 3 keV (A < 4.13 A) is thus to impose some subtle 
downwards curvature on the continuum towards lower energies. 

The only significant absorption features occurring directly in 
the 3 — 10 keV energy range will be due to K-shell transitions in 
the abundant Z > U elements (S, Ar, Ca, Fe and Ni). Absorp- 
tion in the 2 — 5.5 keV range could be due to S, Ar and Ca while 
Fe and Ni may absorb above 6.40 keV and 7.47 keV, respectively. 
Comparison with the high quality Chandra HETGS spectrum of 
the Seyfert 1 galaxy NGC 3783, which shows a very strong warm 
absorber (Kaspi et al. 2002), suggests that absorption due to S, Ar 
and Ca is likely to be weak (see the top panel of their figure 1 and 
also Blustin et al. 2002, Behar et al. 2003). This is confirmed by 
an analysis of the Chandra HETGS spectrum of MCG-6-30-15 (I. 
Lee, piv. comm.; see also section l4'.5.3t . 

Turner et al. (2003a) fitted the RGS spectrum from the 
2001 XMM-Newton observation using a multi-zone warm absorber 
model plus absorption by neutral iron (presumably in the form of 
dust). Although the transmission function of this absorption model 
recovers above a few keV it nevertheless predicts the absorption 
remains significant at higher energies (a ~ 10 per cent effect at 
3keV) due to the combined low-energy edges discussed above. 
Therefore, to account for the absorption in MCG-6-30-15 this 
multi-zone absorption model was included in the spectral fitting, 
with the column densities and ionisation parameters kept fixed at 
the values derived from fits to the RGS data (these are in any case 
very poorly constrained by the EPIC data above 3 keV). The values 
used are tabulated in Table 1 of Turner et al. (2003a; model 2). The 
photoelectric edges were included but not the absorption lines. This 
absorption model therefore accounts for the subtle spectral curva- 
ture imposed by the warm absorber. Resonance absorption lines 
from S, Ar, Ca and Ni are unlikely to contribute significant equiv- 
alent width. The possibility of resonance line absorption by Fe is 
explored separately below. 

4.3.2 Iron K-shell absorption 

The other noteworthy effects of the warm absorbing material above 
3 keV are due to the K-shell of iron. Iron in F-like to H-like ions 
(Fe"'"^'^ — Fe^^^) can absorb through Ka resonance lines in the 
6.4 - 6.9 keV range (Matt 1994; Matt, Fabian & Reynolds 1997) 
which would be poorly resolved (if at all) in the EPIC spectra. Sako 
et al. (2003) estimated the resonance lines could produce a total ab- 
sorption equivalent width of EW ~ -30 eV in MCG-6-30-15 
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based on an analysis of the Fe L-shell transitions present in the 
RGS spectrum from the 2000 XMM-Newton observation. This ab- 
sorption could alter the shape of the observed emission around the 
Fe Kq emission hne and was therefore accounted for in the anal- 
ysis (paper I briefly explored the possible effect of iron resonance 
absorption). 

The other major source of opacity is K-shell photoelectric ab- 
sorption of photons above 7.1 keV (Palmeri et al. 2002). This could 
affect the continuum estimation on the high energy side of the line, 
and therefore the apparent line profile, if not properly accounted 
for (see discussion in e.g. Pounds & Reeves 2002). The edge due to 
neutral iron is negligible in MCG-6-30-15; the Galactic and intrin- 
sic (dust) neutral absorption gave r < 0.01. The warm absorbing 
material contains Fe in wide a range of charge states, leading to 
a blurring together of the many K-edges over the ~ 7. 1 — 9 keV 
range (see Palmeri et al. 2002). In the Turner et al. (2003a) absorp- 
tion model the total optical depth due to these edges (estimated by 
comparing the model transmission below the Fe I edge to that at the 
minimum of the resulting absorption trough) is only Tmax ~ 0.02. 
Using the ionic columns for Fe derived by Sako et al. (2003) gave 
similar results (rmax ~ 0.01 around the K-edges). 

This suggests that while the line-of-sight material does pro- 
duce iron K-edge absorption, it is rather shallow and spread over 
the 7 — 9keV range. The resulting distortion on the transmitted 
spectrum is far smaller than the statistical errors on the data above 
7 keV and thus is of little consequence to the spectral fitting. For 
completeness the absorption due to dust and warm gas in MCG- 
6-30-15 (as derived by Turner et al. 2003a) was included in the 
absorption model used below. Furthermore, the 7—10 keV EPIC 
data were searched for additional Fe absorption (see section U.S.H . 
such as Fe XXV - XXVI edges which would be undetectable in the 
lower energy grating spectra (due to the lack of L-shell electrons). 



4.3.3 Reflection continuum 

Observations of MCG-6-30-15 extending beyond lOkeV with 
Ginga (Pounds et al. 1990; Nandra & Pounds 1994), RXTE (Lee 
et al. 1998, 1999) and BeppoSAX (Guainazzi et al. 1999; paper I) 
revealed an upturn in the spectrum and indicated the presence of 
a reflection continuum. Particularly relevant to the present analy- 
sis is that BeppoSAX observed MCG-6-30-15 simultaneously with 
the 2001 XMM-Newton observation. The data from the PDS instru- 
ment provided spectral information extending up to ~ 100 keV. 
Paper I and Ballantyne et al. (2003) presented the results of fitting 
the BeppoSAX data with reflection models (see in particular section 
3 of Ballantyne et al. 2003). These results strongly indicated that 
the reflection continuum is strong {R ^ 2). 



4.3.4 Narrow iron emission line 

As discussed in SectionQ various missions have resolved the broad 
emission feature spanning the 4 — 7keV range. In each observa- 
tion the spectrum can be fitted in terms of a disc line. There could 
in principle also be a contribution to the reflection spectrum from 
more distant material (S> 100 rg) which would produce a much 
narrower, symmetric line. Such lines have been resolved in the 
spectra of many other Seyfert 1 galaxies (e.g. Yaqoob et al. 2001; 
Yaqoob, George & Turner 2002; Pounds et al. 2003a; Reeves et al. 
2003). The presence of such a narrow 'core' could significantly af- 
fect model fits to the iron line profile if not correctly accounted for 
(e.g. Weaver & Reynolds 1998). 



Table 1. Best-fitting parameters for Gaussian models of the iron line core 
region (5 — 8 keV). ^ indicates the parameter was fixed. 



Line 


£;(keV) 


cr(eV) 


EW(eY) 




1 


6.38 ±0.02 


291 ± 33 


145 ± 15 


318.5/291 


1 
2 


6.41 ±0.03 

6 73+0 01 

"•'■^-O.fM 


312 ± 33 
10-^ 


167 ± 20 

-11 ± 6 


300.2/289 


1 

2 

3 


6.35 ± 0.05 

6 74+0-15 
°-'*-0.09 

6.40/ 


405 ± 50 
10/ 
10/ 


148 ± 20 
^ -8.6 
20 ±5 


263.3/288 



In the case of MCG-6-30-15 the variations in the line pro- 
file observed by ASCA during a low-flux period (the 'deep mini- 
mum') suggested the contribution to the line from distant material 
was small (Iwasawa et al. 1996). The higher resolution Chandra 
HETGS spectrum (Lee et al. 2002) was able to resolve the peak 
of the line emission near 6.4 keV and confirm that the contribu- 
tion due to distant material (i.e. an unresolved line 'core') is weak 
(EW < 20 eV). The 2000 XMM-Newton observation of MCG- 
6-30-15 also showed evidence for narrow line emission (WOl; 
Reynolds et al. 2003). Therefore the possibility of an unresolved 
iron emission line is allowed in the spectral models considered be- 
low. 



4.4 Fitting tlie Fe line core 

Figure |4| shows residuals when the EPIC spectra are fitted with a 
simple model comprising a power-law plus neutral and warm ab- 
sorption (as discussed in section |4.3.1> . The only free parameters 
were the power-law slope (photon index F) and normalisation. The 
spectral interval immediately around the Fe K features (5 — 8 keV) 
was ignored during the fitting, and the model was subsequently 
interpolated over these energies. The residuals clearly reveal the 
strong iron Ka emission line peaking near 6.4 keV. The emission 
line appears asymmetric in this ratio plot, extending down to at 
least ~ 5 keV (this red wing on the line is far broader and stronger 
than that expected from the Compton shoulder from cold gas; Matt 
2002). 

As a first step towards modelling the iron features the 5 — 
8 keV region was examined in detail. A model was built in terms 
of a power-law continuum (absorbed as described in section l4'.3.I> 
plus three Gaussian lines. The best-fitting parameters are shown 
in TableQ The lines were added to model the resolved line emis- 
sion (line I), the unresolved absorption at « 6.7 keV (line 2) and 
the unresolved 6.4 keV emission line (line 3). Modelling the line 
with the single resolved Gaussian emission line provided a reason- 
able fit, but the inclusion of the unresolved absorption and emission 
lines significantly improved the quality of the fit. The inclusion of 
the absorption line was an attempt to account for possible Fe reso- 
nance absorption lines fsection l4.3.2> . 

Figure |5| shows the spectrum around the peak of the Fe emis- 
sion fitted with the three line model (only EPIC pn:s data are shown 
but all four spectra were included in the fitting). The resolved emis- 
sion line accounted for the majority of the flux from the core of 
the Fe emission line, and was significantly resolved with a veloc- 
ity width FWHM « 4.5 x 10* kms"^ . Fitting the four EPIC 
spectra individually with the same model gave consistent results; in 
each case the emission line was very significantly resolved and the 
measured widths were consistent with one another. 
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Figure 4. Ratios of EPIC spectra to absorbed power-law model. The model 
was fitted to the 3 — 5 and 8 — 10 kcV data. Note the model underestimates 
the flux on the red side of the line (~ 5 — 6.4 koV) while the flux is overes- 
timated above 7 keV. Also note the narrow 'notch' at ^ 6.7 keV noticeable 
in the pn:s spectrum. 



The high dispersion velocity width of the emission line core 
is robust to the details of the underlying continuum fitted over this 
fairly narrow energy range (it remain significantly resolved after 
including additional cold gas and/or Fe K edge absorption in the 
model; see also section lT.S.H . Forcing the entire emission line to be 
narrow (i.e. unresolved) gave an unacceptable fit (Ax^ = +266.2 
compared with the resolved line). In addition to ~ 6.4 keV emis- 
sion there could also be Fe lines at energies ~ 6.7 and ~ 6.9 keV, 
corresponding He-like and H-like ions, respectively. The simul- 
taneous presence of iron lines at these three energies (convolved 
through the EPIC response) could in principle conspire to make the 
Fe emission mimic a single, resolved line. See Matt et al. (2001) 
and Bianchi et al. (2003) for an examples where such line blends 
have been observed in EPIC spectra of Seyferts. Allowing for the 
presence of three unresolved lines (at energies of ~ 6.4, ~ 6.7 
and ~ 6.9 keV) also gave an unacceptable fit. The Fe emission was 
thus unambiguously resolved by EPIC. 

This is of great interest because the resolved line core is con- 
siderably broader than those of some other Seyfert I galaxies (e.g. 
Yaqoob et al. 2001; Kaspi et al. 2002; Pounds & Reeves 2002; 
Page et al. 2003a; Page et al. 2003c) and suggests an origin close 
to the central SMBH. Neglecting for the time being the contribu- 
tion of any strongly redshifted component to the line profile, the 
width of the resolved core can be used to infer the radial distance 
of the line emitting gas under the assumption that the material is 
gravitationally bound (and emits at 6.40 keV in the rest frame). 



5 5.5 6 6.5 7 7.5 B 

Energy (keV) 

Figure 5. EPIC pn:s spectrum immediately around the core of the Fe line 
(crosses). The histograms show the spectral model components, namely a 
power-law continuum, a broad (cr fa 410 eV) Gaussian emission line at 
« 6.4 keV. an unresolved emission line at 6.4 keV and an unresolved ab- 
sorption feature at 6.7keV. 



In this case the rms velocity width of the line (a) is determined 
by its proximity to the SMBH (r/rg) and the geometry according 
to: r/rg ^ {c/a)^/q (where q depends on the geometry; Krolik 
2001). Assuming randomly oriented circular orbits (q — 3) places 
the line emitting material at r ~ 80 rg, while assuming a Keplerian 
disc (inclined at i = 30°; q = 8) requires the material to be at 
r ~ 30 Tg. The resolved line core thus suggests emission from near 
the central SMBH, necessitating the use of models accounting for 
the relativistic effects operating in this region. 

Fitting this model over the 3 — 10 keV range gave a bad fit 
(X^ = 983.1/777 dof) with the emission line energy fixed to 
6.4 keV. Allowing the energy of the broad Gaussian to be free 
improved the fit (x^ = 862.1/776 dof) but the Gaussian be- 
came extremely broad and redshifted (E — 5.9 ± 0.2 keV and 
a ^ 0.7keV). Fitting with two Gaussians to model the broad 
line (to account for the resolved core and also a highly red- 
shifted component), as well as the narrow 6.4 keV emission and 
« 6.7 keV absorption, provided a good fit (x^ = 730.7/773 dof). 
The two broad Gaussians had the following parameters: Ei — 
6.38tJ! °| keV, cri = 352j:^°^eV, EWi ^ 120 eV and E2 = 
4.9l?;^'keV, (72 = 1.0l?,:^keV, EW2 ^ 150 eV. The energy 
and width of the second Gaussian further suggests the presence of 
a highly broadened and redshifted component to the Fe emission. 
However, such a claim is dependent on the assumed form of the 
underlying continuum. 



4.5 Dependence on the continuum model 

The simple absorbed power-law continuum model used above 
(Fig-EJ clearly overestimated the continuum above 7keV (on the 
'blue' side of the line). These residuals are a direct result of fit- 
ting the continuum with an over-simplistic model, and specifically 
could be caused by of one of the following: (a) additional K-shell 
absorption by Fe, (6) intrinsic curvature in the continuum, (c) a 
continuum curved by excess absorption, or (d) excess emission on 
the red side of the line (3 — 5.5 keV). These last two possibilities 
would lead to the true continuum slope over the 3 — 10 keV range 
being underestimated (the fit was driven by the better statistics in 
the 3 — 5 keV region). 

The first three of these alternative possibilities (a — c), illus- 
trated in figure|6| are briefly discussed below and found to be unsat- 
isfactory. In the subsequent analysis of sections l4!6l and l477l it was 
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Figure 6. Schematic of three alternative models that do not include an 
emission feature on the red side of the iron Hne. The models comprise a 
power-law continuum with resolved Gaussian line in addition to: (a) addi- 
tional Fe K-edge absorption, (b) a break in the power-law slope at ss 6 keV 
and (c) enhanced soft X-ray absorption. 
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Figure 7. Confidence contours (68.3, 90, 99 per cent limits) on absorption 
edge parameters when the 7—10 kcV spectrum was searched for absorp- 
tion. 



assumed that there is excess emission extending from the Fe line 
to lower energies (i.e. possibility d), and this was modelled using 
relativistic and non-relativistic emission features, respectively. 



4. 5. 1 Additional iron K absorption ? 

As discussed in section 14.3.21 the Fe K-shell absorption edges ex- 
pected based on the grating spectra has been included in the ab- 
sorption model. It remains possible that there is absorption from 
additional Fe that was not accounted for in this model. 

The 7 — lOkeV EPIC spectrum was searched for additional 
absorption by fitting the data over only this limited energy range 
with a power-law and an absorption edge. Although this model is 
very simple, over the limited energy range the underlying contin- 
uum should not deviate noticeably from a power-law, and with the 
limited resolution and statistics available at these energies an edge 
can provide an approximate description of additional Fe absorption. 

The data constrain the depth of any such edge to be r < 0.06 
(90 per cent confidence) over the energies expected for Fe K 
(7.11 — 9.28 keV). In particular the depth of edge allowed at the 



energies expected for Fe XXV - XXVI (at the rest-frame velocity 
of MCG-6-30-15) is r <J 0.03, ruling out the presence of a large 
column of such gas (which would have gone undetected in the low 
energy grating spectra). 

The best fitting edge model gave E — 7.31 ± 0.14keV and 
r = 0.036 ±0.024. Figure|7|shows the confidence contours for the 
edge parameters in this 'sliding edge' search. Including this edge 
in the model made virtually no difference to the residuals around 
the iron line. In particular the apparent asymmetry remained. That 
said, the best-fitting optical depth should perhaps be treated as an 
upper limit due to the likely presence of Fe K/3 emission at i5 ~ 
7.06 keV, which must accompany the Kq emission. Allowing for 
the possibility of line emission al E = 7.06 keV resulted in the 
edge depth becoming consistent with zero, with a limit of r ^ 
0.06. Including the Fe K edge absorption contained in the model 
discussed in section |4.3 . 1 1 further reduced the depth of any excess 
Fe absorption. 

A caveat is that if a large column of Fe is present with a range 
of intermediate charge states the resulting absorption will not re- 
semble an edge (Palmeri et al. 2002). However, the low energy 
grating data do not support the presence of such a large column 
of ionised Fe in addition to that already included in the model (Lee 
et al. 2001; Sako et al. 2003; Turner et al. 2003a). 



4.5.2 An intrinsically curved 3 — 10 keV continuum? 

Ballantyne et al. (2003) and paper I presented analyses of the simul- 
taneous BeppoSAX PDS spectrum of MCG-6-30-15. After allow- 
ing for the reflection continuum, the underlying power-law did not 
show any evidence for intrinsic curvature. In particular, the lower 
limit on the energy of an exponential cut-off in the continuum was 
constrained to be Ecut ^ 100 keV. Thus any subtle spectral cur- 
vature due to the continuum deviating from a power-law at high 
energies (as predicted by thermal Comptonisation models) will not 
affect the 3—10 keV bandpass. 

As a check for the possible effect of more drastic curvature in 
the 3 — 10 keV band, the EPIC spectrum was fitted with a broken 
power-law model, after excluding the region immediately around 
the Fe K features (5-8 keV). This provided a good fit (xl = 0.92). 

'i.Sbtlil keV where 
above 



The break in the continuum occurred at i? 
the slope changed from Ti = 2.00^o 04 to r2 = 2^14+° °^ 
the break energy (the slopes are given based on the pn:s spectrum, 
similar differences were derived from the other spectra). This form 
of spectral break, with the spectral slope increasing by AF ~ 0.15 
above 5 keV is highly unusual and unexpected when compared to 
other Seyfert I spectra. Furthermore, the strong reflection compo- 
nent seen in the BeppoSAX PDS data (section l4.3.3> should lead to 
a flattening of the spectral slope at higher energies. Indeed, when 
the broken power-law model was interpolated across the 5 — 8 keV 
region it predicted a broad iron line extending down to at least 
5 keV (the energy of the continuum break) but appeared to match 
poorly the continuum slope at energies above the line (the model 
was too steep compared to the data above 7keV). An intrinsically 
curved/broken 3—10 keV continuum is therefore considered highly 
unlikely. 



4.5.3 Excess soft X-ray absorption ? 

The warm absorbing gas could perhaps impose curvature on the 
3 — 5.5 keV continuum in excess of that already accounted for in 
the absorption model (section l4.3.1> . Significantly higher columns 
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Figure 8. Ratio of EPIC pn:s spectrum to an absorbed power-law model. 
The model was fitted over the ranges 3 — 5keV and 8 — lOkeV with the 
absolution model discussed in section l4.3.1l but allowing the parameters 
of one of the absorbing zones to be free (see section l4.5.3l . There is little 
difference in the pattern of residuals compared with figure |4] In particular 
there is excess emission on the red side of the line extending down to at least 
5keV, and again the model slightly overestimates the flux above 7koV. 



of O or Si ions, for example, would have the effect of increasing 
the opacity in the high energy tail of the absorption edges occur- 
ring at energies <^ 3 keV. This can be quite effectively modelled 
by simply allowing the parameters of one of the warm absorbing 
zones (namely the column density and ionisation parameter) to be 
free in the fitting. The data were therefore fitted in the 3—10 keV 
range, again excluding Fe K region (5 — 8 keV), using a power-law 
modified by the absorption model of section l4.3.1l but allowing the 
parameters of one of the warm absorber zones to vary. The result 
was that column density of the low ionisation absorber increased to 
fit better the slight curvature in the 3 — 5 ke V band. The fit was good 
(Xv — 0.92) but when the model was interpolated into the Fe K re- 
gion the residuals again suggested a strongly asymmetric emission 
feature (as shown in Figure |8}. The inclusion of extra absorption 
below 3 keV did therefore not alter the requirement for an asym- 
metric, red emission feature. However, this model is disfavoured 
as it severely under-predicted the flux level at lower energies when 
extrapolated below 3 keV. Very similar results were obtained af- 
ter allowing the other (higher ionisation) absorbing zones to vary 
during the fitting. 

The other possibility, namely that the 3 — 5.5 keV spectrum 
is significantly affected by S, Ar and/or Ca absorption, was also 
explored. Absorption by the He- and H-like states of these ions 
will not affect the spectrum < 2 keV (except if they are accom- 
panied by significant columns of M- and/or L-shell ions; Behar & 
Netzer 2002) and therefore the soft X-ray EPIC and RGS data do 
not strongly constrain their presence. The upper limit on the opti- 
cal depth of an absorption edge is Tmax < 0.02 throughout most 
of this energy range. The only exception is at ~ 4 keV where the 
addition of a weak edge (r ~ 0.03) did significantly improve the 
fit. Thus it remains plausible that some absorption by, for exam- 
ple, highly ionised Argon (e.g. Ar XVII) may weakly contribute in 
this region. However, this feature must be considered with caution. 
An examination of the EPIC calibration data for 3C 273 (see Ap- 
pendix]^ revealed similar residuals at the same (observed frame) 
energy. This feature could thus be due to an instrumental effect at 
the « 3 per cent level (possibly due to Ca residuals on the mir- 
ror surfaces; F. Haberl priv. comm.). In any event the affect on the 
derived iron line parameters was negligible. 



4.6 Models including strong gravity 

The spectral fits detailed above demonstrated that the Fe emission 
line was significantly resolved, implying emission from near the 
SMBH. These did not directly address the issue of a highly red- 
shifted component to the line. Furthermore, MCG-6-30-15 shows 
a strong reflection continuum (section l4.3.3t which needs to be ac- 
counted for in the spectral modelling. Therefore, in this section the 
3 — 10 keV spectrum is fitted with models accounting for emis- 
sion from an X-ray illuminated accretion disc extending close to 
the SMBH (see also WOl, paper I and references therein). 

The code of Ross & Fabian (1993; see also Ballantyne et al. 
2001) was used to compute the spectrum of emission from a pho- 
toionised accretion disc. The input parameters for the model were 
as follows: the photon index (F) and normalisation (TV) of the in- 
cident continuum, the relative strength of the reflection compared 
to directly observed continuum (R) and the ionisation parameter 
(^ = An Fx /nil)- This model computes the emission (including 
relevant Fe Ka lines), absorption and Compton scattering expected 
from the disc in thermal and ionisation equilibrium. 

To account for Doppler and gravitational effects, the disc 
emission spectrum was convolved with a relativistic kernel cal- 
culated using the Laor (1991) model. The calculation was carried 
out in the Kerr metric appropriate for the case of a disc about a 
maximally spinning ('Kerr') black hole. The model uses a broken 
power-law to approximate the radial emissivity (as used in paper I). 
The parameters of the kernel were as follows: inner and outer radii 
of the disc (rin, rout), inclination angle (i), and three parameters 
describing the emissivity profile of the disc, namely the indices 
(<Zin,gout) inside and outside the break radius (rbr)- In addition, 
narrow Gaussians were included at « 6.4 keV to model a weak, 
unresolved emission line and at ~ 6.7 keV to model the unresolved 
absorption feature. This emission spectrum (power-law, reflection, 
narrow emission line) was absorbed using the model discussed in 
section l4.3.1l 

This model gave an excellent fit to the data (x^ ~ 
652.3/771 dof). The residuals are shown in Fig.|9| The best-fitting 
parameters corresponded to strong reflection (R = 1.48lQ Qg) 
from a weakly ionised disc (log(^) <. 1.4). The parameters of the 
blurring kernel defined a relativistic profile from a disc extending 
in to Tin — 1.8 ± 0.1 Tg with a steep (i.e. very centrally focused) 
emissivity function within « 3rg and a flatter emissivity with- 
out. Specifically, the emissivity parameters were gin — 6.9 ± 0.6, 
Qout = 3.0 ± 0.1 and rbr = 3.4 ± 0.2 rg and the inclination of 
the disc was i — 33 ± 1°. The narrow emission and absorption 
lines at 6.4 keV and « 6.7 keV had equivalent widths of ~ 10 eV 
and « — 10 eV, respectively. This model is very similar to the best- 
fitting disc model obtained in paper I from the MOS and BeppoSAX 
data. The only notable difference between the two models is that the 
emissivity is more centrally focused in the present model. This is 
primarily a consequence of including of the EPIC pn data in the fit- 
ting, which better defined the red wing of the line, and allowing for 
narrow Fe Kq absorption and emission lines, which subtly altered 
the appearance of the line profile. 

The above model includes emission from well within 
6rg, the innermost stable, circular orbit (ISCO) for a non- 
rotating ('Schwarzschild') black hole. In order to test whether a 
Schwarzschild black hole (or specifically rin ^ 6rg) is incompat- 
ible with the data the ionised disc model was refitted after convo- 
lution with a diskline kernel appropriate for emission around 
such a black hole (Fabian et al. 1989). The disc was assumed to 
be weakly ionised (log(^) — 1), the inner radius was fixed at 
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Figure 9. Residuals from fitting the EPIC spectra with a variety of models 
(for clarity only the EPIC pn:s data are shown). The models comprise an 
absorbed power-law continuum plus: (a) broad Gaussian emission line, (b) 
two broad Gaussians, (c) relativistically blun'ed emission from an accretion 
disc and (d) broad Gaussian and partially covered continuum. Also included 
is all these models are an unresolved absorption feature at fa 6.7keV and 
an unresolved, neutral iron emission line at 6.4 keV. 



^in = 6 rg in the fitting, and the emissivity was described by an un- 
broken power-law. This model provided an acceptable fit to the data 
(X^ = 698.9/775 dof) albeit substantially worse than the model 
including emission from within 6 rg (Ax^ = 46.6 between the two 
models). Thus a model in which the broad iron line and reflection 
spectrum originate from an accretion disc about a Schwarzschild 
black hole cannot be ruled out based on the 3 — 10 keV EPIC spec- 
trum. However, the model including emission down to rin ~ 1-8 rg 
was preferred for two reasons. Firstly it provided a much better fit 
to the data. Secondly the derived strength of the reflection spectrum 
for the Tin ~ 1.8 rg model (R ~ 1.4 — 1.8) was more compatible 
with that derived from the independent analysis of the BeppoSAX 
PDS data (R ^ 2; section I4.3.3> . The best-fitting value for the 
Schwarzschild model was only R — 0.77 ± 0.04. 



4.7 Models excluding strong gravity 

The previous section discussed models in which the spectrum is 
shaped by strong gravity effects. In this section some trial spectral 
models are discussed that attempt to explain the spectrum without 
recourse to strong gravity. In particular, the core of the iron line 
emission is modelled in terms of a resolved Gaussian and an un- 
resolved emission line, both centred on 6.4 keV (see section |4^ . 



This approximates the distortion on the Fe emission line due to 
Doppler but not gravitational effects. The asymmetry around the 
line, in particular the broad excess on the red side, is instead mod- 
elled in terms of additional emission components rather than al- 
lowing for highly redshifted iron line emission. The possibilities 
explored are as follows: (i) a blackbody, (ii) a blend of broadened 
emission lines from elements other than Fe, and (iii) a partially 
covering absorber. 

(i) Blackbody emission could provide a broad bump not dissimi- 
lar from an extremely broad line. A model comprising a power-law, 
three Gaussians (resolved and unresolved Fe emission lines and the 
unresolved absorption feature at ~ 6.7 keV) and a blackbody gave 
a reasonable fit to the data (x^ = 767.0/775 dof). The blackbody 
temperature was kT « 1.1 keV and it produced a broad excess in 
flux from < 2 keV to ~ 6 keV, i.e. on the red side of the line. This 
model is therefore formally acceptable but the fit is considerably 
worse than with the relativistic disc models. 

(ii) Along with the Kq emission from Fe there may also be Ka 
emission from e.g. S, Ar, Ca or Cr that would also be Doppler 
broadened. An additional broad Gaussian line was included in the 
model with a centroid energy in the range 3 — 5.5 keV, and a width 
fixed to be the same as the resolved Fe line. This provided an unac- 
ceptable fit (x^ = 920.6/780 dof). Adding a third broad Gaussian 
yielded an acceptable fit (x^ = 809.7/778 dof). The additional 
lines were centred on ~ 3.7 keV and ~ 5.0 keV. The first of these 
could plausibly be due to Ca XIX-XX, but may be caused by a cali- 
bration artifact (see section U.S.Bt . No strong X-ray line is expected 
at ~ 5.0 keV. The line equivalent widths were ~ 105 keV and 
~ 60 keV, at least an order of magnitude higher than expected for 
fluorescence from low-Z elements (e.g. Matt et al. 1997; Ross & 
Fabian 1993). 

(iii) The presence of partially covered emission (i.e. a patchy 
absorber; Holt et al. 1980) can cause 'humps' in the 2 — lOkeV 
spectral range. In such model the observed spectrum is the sum 
of absorbed and unabsorbed emission. A neutral partial covering 
absorber (with column density, A^h, and covering fraction, /c, left 
as free parameters) was applied to the power-law continuum in the 
spectral model. Such a model can provide a reasonable fit to the 
data (x^ = 745.7/775 dof) with Nn ^ 1.4 x 10^^ cm"^ and 
/c = 0.27 ±0.03. 

This analysis demonstrated that alternative emission compo- 
nents can reproduce the broad spectral feature observed in the 
3 — lOkeV EPIC spectrum. However, these models all give best- 
fitting x^ values considerably worse than the relativistic disc model 
discussed in section 1461 Possibilities (i) and (ii) are rather ad hoc 
and physically implausible. A partially covering absorber seems 
somewhat more reasonable since partial covering has been claimed 
for other Seyfert galaxies. However, such a model cannot explain 
the high energy data, in particular, the need for a strong reflection 
component (section l4.3. Bl and see also Reynolds et al. 2003). Fur- 
thermore, allowing for partial covering in the model did not elimi- 
nate the need for a strongly Doppler broadened iron line. 

As a final test of the possible impact of partial covering on 
fitting the relativistic disc models, the model from section l46l was 
re-fitted including partial covering. The addition of a partially cov- 
ering absorber to the best-fitting reflection model did not provide a 
significant improvement to the fit (Ax^ ~ 2.6 for two additional 
free parameters). The derived covering fraction of the absorber was 
poorly constrained (0.21 ^ /c ^ 1.0) with a column density 
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Figure 10. Top panel: correlation diagram comparing EPIC pn count rates 
in two different energy bands. The points are from 1000 s binned light 
curves. The two observations (2000 and 2001) are marked with different 
symbols. The dashed hne represents the best fitting linear function. Bottom 
panel: correlation diagram after averaging the data in bins of 20 points. The 
points are shown with error bars on both axes, which are of comparable size 
to the points. The dashed line represents the best fit to the binned data. 



profile parameters (specifically rin = 1.8 ± 0.1 rg). This analysis 
suggests that the possible presence of partially covering absorption 
does not alter the derived parameters for the relativistic emission 
line. Furthermore, the inclusion of the partial coverer resulted in a 
lower value for the reflection strength (R ~ 1.1), in contradiction 
with the PDS spectrum (ii 2 see section U.B.St . 



5 SPECTRAL VARIABILITY 

The rapid X-ray variability of MCG-6-30-15 has long been known 
to show complex energy dependence (e.g. Nandra, Pounds & Stew- 
art 1990; Matsuoka et al. 1990; Lee et al. 1999; Vaughan & Edel- 
son 2001). Vaughan et al. (2003a) used Fourier methods to examine 
differences between broad energy bands as a function of timescale. 
In this section the spectral variability is examined in a finer en- 
ergy scale, at the expense of cruder timescale resolution, in order 
to probe the details of changes due to specific spectral components. 
The EPIC data were used over the 0.2—10 keV band as the calibra- 
tion issues discussed above (see also Appendix IaI will not affect 
these analyses. 

5.1 Flux-flux analysis 

In this section the spectral variability is investigated using 'flux- 
flux' plots. Taylor, Uttley & M'^Hardy (2003) used the RXTE data 
of MCG-6-30-15 to investigate the correlations between the fluxes 
in different energy ranges. This section parallels their analysis. 

Figure [Tol shows the count rates in a soft band (1 — 2keV) 
against the simultaneous count rates in a hard band (3 — 10 keV) 
from the 1000 s binned EPIC pn light curves. Clearly the two sets 
of count rates (hereafter fluxes) show a very significant correlation. 
This is true for both XMM-Newton observations. The data were fit- 
ted with a linear model (of the form F{E2) = mF{Ei) + c; where 
F{Ei) and F{E2) represent the fluxes in the two bands) account- 
ing for errors in both axes. This resulted in a formally unacceptable 
fit (rejection probability > 99.99 per cent). However, the poorness 



of the fit is a result of intrinsic scatter in the correlation, rather than 
any non-linearity. This intrinsic scatter in the relation between dif- 
ferent energy ranges is a manifestation of both the difference in 
shape of the power spectral density (PSD) function and low coher- 
ence between variations in different energy bands (Vaughan et al. 
2003a). Therefore, to overcome this the data were then binned such 
that each bin contains 20 points and errors were calculated in the 
usual fashion (equation 4.14 of Bevington & Robinson 1992). The 
binned data, which reveal the average flux-flux relation, provided a 
reasonable fit to a linear function (rejection probabilities of 33 and 
89 per cent for the 2000 and 2001 datasets, respectively). 

As concluded by Taylor et al. (2003) the tight linearity of the 
flux-flux correlation (at least over the relatively short time scales 
probed by each XMM-Newton observation) strongly indicates that 
the flux variations are dominated by changes in the normalisation 
of a spectral component with constant softness ratio m (i.e. the gra- 
dient of the flux-flux correlation). In addition there must be an ad- 
ditional spectral component that varies little and contributes more 
in the 3 — 10 keV band than the 1 — 2 keV band and therefore pro- 
duces the positive constant offset c. In other words the flux in each 
band is the sum of a variable component and a constant component: 



F{E) = NvV{E) + NcC{E) 



(1) 



(where AV and A^c are the normalisations of the variable and con- 
stant spectra V{E) and C{E) respectively) ^. The fact that the re- 
lation between F{Ei) and F{E2) is linear means that the gradient 
gives m = V{E2)/V{Ei) and the offset gives c = C{E2) - 
mC{E\). As can be seen from the figure the relation changed be- 
tween the 2000 and the 2001 observations; these changes in slope 
and offset imply the shape of the two spectral components dif- 
fered slightly between the two XMM-Newton observations. The 
flux-correlated changes in the spectrum are caused by changes in 
the relative contributions of these two spectral components. This 
'two component' interpretation of the X-ray spectral variability of 
Seyfert Is has been discussed by e.g. Matsuoka et al. (1990); Nan- 
dra (1991); M'Hardy et al. (1998); Shih et al. (2002); Fabian & 
Vaughan (2003). 



5.1.1 Decomposing the spectrum 

As discussed in Fabian & Vaughan (2003), the 'difference spec- 
trum' can be used to determine the shape of the variable component 
of the spectrum, V{E) (see section l5!2l . The offset of the flux-flux 
plot, c, gives a measure of the strength of the constant component in 
the spectrum and so by measuring the offset as a function of energy, 
c{E) it is possible to estimate the spectrum of the constant compo- 
nent C{E) (see section 4.3 of Taylor et al. 2003). Figure fTTI shows 
the spectrum of the constant component in MCG-6-30-15 deduced 
by measuring the offset c{E) of the flux-flux relations. Binned flux- 
flux plots were constructed by comparing the light curves in each 
energy band to the 1.0 — 1.5 keV band {F{E-i)). In each case the 
constant offset was measured from the best fitting linear model and 
normalised to the average flux in the energy band under examina- 
tion. This yielded an estimate of the fractional contribution to the 
spectrum from non-varying components. 

The spectrum of the constant component shown in figure fTTI 
reveals that the constant component is soft below ~ 1 keV and 
hard at higher energies. The largest fractional contribution from the 



The spectra of the two components V{E) and C{E) include the effects 
of absorption. 
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Figure 11. Fractional contribution to tlie spectrum of tlie constant compo- 
nent {C{E)/ F{E) deduced from the linear flux-flux relation). The en'ors 
are la confidence limits obtained from the linear fit to the flux-flux data. If 
the constant component contributes a non-zero fraction in the 1.0 — 1.5 keV 
range the contributions at all energies will increase but the spectral shape 
will remain the same. The arrows indicate the absolute upper limits on the 
fraction of constant emission (see text). 



constant component occurs in the iron K band, where it contributes 
~ 40 per cent of the total flux. There are two assumptions im- 
plicit in this analysis. The first is that the variable component has a 
constant spectral shape V{E) when averaged as a function of flux 
(which is implied by the linearity of the flux-flux relations). The 
second assumption is that there is a negligible constant component 
in the 1.0 — 1.5 keV band (which was used as the comparison), i.e. 
C{Ei) = 0. If the 1.0 — 1.5 keV emission does contain a non- 
zero constant component (C{Ei) > 0) then the absolute scale of 
the inferred spectrum would increase but its spectral form will be 
approximately the same. The maximum possible amount of con- 
stant emission in the 1.0 — 1.5 keV band is given by the minimum 
flux in that band (i.e. C{Ei) < mm[F{Ei)]), therefore this 
can be used to place an absolute limit on the strength of the con- 
stant component: c{E) ^ C{E) ^ {c{E) + m{E) mm[F{Ei)]}. 
This upper limit is also indicated in figure [TTI It should be noted 
that the constant component C{E) must be affected by the soft X- 
ray absorption. The lack of features associated with absorption in 
figure fm fin particular the lack of a spectral jump at ~ 0.7 keV) 
demonstrates that the constant component is affected by absorption 
exactly as is the total spectrum (hence the absorption was factored 
out by normalising the spectrum C{E) by the total spectrum). This 
spectrum can be identified with the Reflection-Dominated Compo- 
nent (RDC) of Fabian & Vaughan (2003) and is henceforth referred 
to as the RDC. 



5.2 Flux-resolved spectra 

A more conventional way to examine how the spectrum of the 
source evolves with flux is to extract energy spectra from different 
flux 'slices.' This was done for the 2001 XMM-Newton data by di- 
viding the pn data into ten flux intervals such that the total number 
of photons collected from the source in each flux interval was com- 
parable (~ 5 X 10^). The highest and lowest flux slices differ by a 
factor of 3 in count rate and all ten are marked on the light curve 
in Figure [T2I The ten spectra extracted from each flux slice were 
then binned in an identical fashion. There were clear, systematic 
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Figure 12. 0.2 - 10 kcV EPIC pn light curve of the 2001 observation. The 
light curve has been divided into ten count rate intervals marked /I — /lO. 

spectral changes between the ten flux intervals, these are illustrated 
in two different ways in figures f73l and ll4l Figure^|shows the ten 
spectra as ratios to a (F = 2) power-law model (modified only by 
Galactic absorption). This illustrates the softening of the 2 — 10 keV 
spectrum as the source gets brighter. In addition the strength of the 
iron line relative to the continuum gets weaker as the source bright- 
ens (i.e. the equivalent width decreases with increasing continuum 
luminosity) meaning that the iron line becomes a lower contrast 
feature at high fluxes. 

Figure^|shows nine of the ten spectra as a ratio to the high- 
est flux spectrum (/lO). This again shows the spectrum became 
harder (in the 2 — lOkeV range) and the iron line became rela- 
tively stronger as the source became fainter. It is also interesting to 
note that at energies below ^ 2 keV the spectrum becomes softer 
as the source becomes fainter, contrary to the trend at higher ener- 
gies. In addition, there is no strong feature in the spectral ratios at 
~ 0.7 keV, indicating that the fractional strength of the ~ 0.7 keV 
spectral jump remains constant with flux. This is consistent with 
the jump being due to absorption with a constant optical depth. 

These average changes in the detailed spectral shape as a 
function of flux can be explained using the simple two-component 
model discussed above. The inferred constant emission comprises 
a soft excess below ~ 1 keV, a harder tail above ^ 2 keV, and 
a strong iron line. One would expect the total (variable+constant) 
spectrum to change exactly as observed: as the total flux decreases 
the contribution from the constant component, relative to the vari- 
able component, increases. Therefore, at low fluxes, the spectrum 
will become softer below ^ 1 keV, harder above ~ 2 keV, and 
display a more prominent iron line. 



5.2.1 Average EPIC difference spectrum 

The five highest flux spectra were combined, as were the five lowest 
flux spectra, to produce two spectra representing the source during 
flux intervals above and below the mean (hereafter 'high' and 'low' 
fluxes). Subtracting the low flux spectrum from the high flux spec- 
trum will remove the contribution from the constant component. As 
discussed in Fabian & Vaughan (2003), under the assumption that 
the total spectrum can be accurately described (on average) by the 
superposition of two emission components, one variable and one 
not (both modified by absorption), the difference spectrum will give 
the spectrum of only the variable component (modified by absorp- 
tion): Fhi{E) - F^E) = (iVy.i^i - iVv.io)l^(£) (where V{E) 
includes absorption). 

Figure [Tsl shows the EPIC pn difference spectrum (compared 
to a power-law model fitted in the 3 — 10 keV range) from the 2001 
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Figure 13. EPIC pn spectra from each of the count rate shces (shown in 
figure IT2I compared to a F = 2 power-law model (modified by Galactic 
absoiption). The data have been shifted upwards to aid clarity. It can be 
seen that as the flux decreases from /lO to /I the overall spectrum becomes 
harder (above 2 keV) and the strength of the iron feature decreases (relative 
to the continuum model). 



observation. A power-law (modified only by Galactic absorption) 
gave an excellent fit to the 3 - lOkeV data (x^ = 28.1/60 dof) 
and yielded a photon index of F = 2.20 ± 0.05. Extrapolating 
this power-law down to lower energies reveals the signature of 
strong absorption, although the exact amount depends on the as- 
sumed shape of the underlying continuum (as illustrated by the up- 
per and lower datasets shown in the figure). It also be noted that if 
the variable spectral component steepens at low energies (i.e. con- 
tains a soft excess) the inferred absorption profile at lower energies 
could be systematically underestimated even further (see section 6 
of Turner et al. 2003a). Leaving aside these uncertainties on the in- 
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Figure 14. EPIC pn spectra from the count rate slices shown as a ratio to 
the highest count rate spectrum (/lO). Similar to figure fTsl the spectra can 
be seen to become harder and show relatively stronger iron features as the 
count rate decreases. 



ferred spectrum at lower energies, it remains true that above 3 keV 
the average spectrum of the varying component V{E) can be de- 
scribed accurately as a single power-law (T « 2.2). In particular it 
should be noted that the difference spectrum (unlike the spectrum 
of the constant component described above) does not possess any 
obvious features in the iron K region. The upper limit on the flux 
of the line core (section|^} present in the difference spectrum is a 
factor ~ 20 lower than that present in the time averaged spectrum. 

A further investigation of the flux-resolved spectra revealed 
that the difference spectrum is consistent with a power-law at all 
flux levels. Nine difference spectra were calculated by subtracting 
the lowest flux spectrum (/I) from each of the nine higher flux 
spectra (/2 — /lO). These nine difference spectra were fitted over 
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Figure 15. EPIC pn difference spectra produced by subtracting the low flux 
spectrum from the high flux spectrum. The difference spectrum is shown as 
a ratio to a power-law (modified by Galactic absoiption) fitted across the 
3 — lOkeV range. The data/model ratio calculated using the best fitting 
power-law slope is shown in black (circles). Above and below this lie the 
ratios calculated assuming the 90 per cent lower and upper Umits on the 
slope of the 3 — 10 keV power-law. The data were rebinned for clarity. 



the 3—10 keV range with an absorbed power-law model. In all nine 
cases the simple power-law model was found to provide a good fit 
to the data (xt ~ 1-0), with no evidence for strong systematic 
residuals around the Fe K band. The limit on the depth of possi- 
ble Fe edges in the 7.1 — 8keV region was r < 0.07, consistent 
with the limits obtained in section 14.5.11 In the discussion below 
this variable component to the spectrum will be referred to as the 
variable Power-law Component (PLC; Fabian & Vaughan 2003). 



5.3 rms spectra 

The root mean squared (rms) spectrum measures the variability am- 
plitude as a function of energy and can in principle reveal which 
spectral components are associated with the strongest variability 
(see Edelson et al. 2002 and Vaughan et al. 2003b). Figure fT6l shows 
two rms spectra calculated from the 2001 XMM-Newton observa- 
tion of MCG-6-30-15, clearly showing in both cases a strong de- 
pendence of the variability amplitude on energy. 

The two rms spectra show the variability amplitude on differ- 
ent timescales. The upper spectrum shows the fractional rms ampli- 
tude integrated over the entire observation in the standard fashion 
(eqn. 1 of Edelson et al. 2002) using 1000 s binned light curves. 
As the variability is dominated by long timescale changes (Uttley 
et al. 2002; Vaughan et al. 2003a) this therefore reveals the en- 
ergy dependence of variations occurring on timescales comparable 
to the length of the observation (~ 100 ks). The lower spectrum 
shows the fractional rms of the point-to-point deviation (i.e. the 
rms difference between adjacent time bins, as defined by eqn. 3 
of Edelson et al. 2002). This therefore only measures fluctuations 
between neighbouring time bins and so probes the energy depen- 
dence of the variability on short timescales comparable to the bin 
size (~ 1 ks). The errors were calculated as in eqn. 2 of Edelson et 
al. (2002), but as discussed in their appendix should be considered 
only as approximations since they strictly assume the light curves 
are drawn from independent Gaussian processes. Matsumoto et al. 
(2003) previously used the longest ASCA observation of MCG~6- 
30-15 to obtain rms spectra on different timescales. 
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Figure 16. Top panel: fractional rms variability amplitudes as a function 
of energy on two different timescales. The upper data represent the ampli- 
tude integrated over the entire (2001) observation while the lower data rep- 
resent the rms deviation between neighbouring points (point-to-point rms) 
and therefore samples only short timescale variability (~ 1000 s). The er- 
rors were calculated as described in the text. The dashed line indicates the 
energy of the (neutral) Fe Ka line (6.40 keV). The histogram shows the 
'two component' model rms spectrum. Bottom panel: The ratio of the two 
rms spectra, showing how the energy dependence of the variability ampli- 
tude changes with variability timescale. 



As is clear from the figures the rms spectra show broad peak 
at 1 keV and a localised depression around the energy of the 
iron K emission line, strongly suggesting a suppression of the vari- 
ability due to the presence of the emission line (see also Inoue & 
Matsumoto 2001). The exact energy dependence of the variability 
amplitude differs between the two timescales (evident from the ra- 
tio of the two rms spectra), as should be expected if the shape of 
the PSD is a function of energy (Vaughan et al. 2003a). 



5.3.1 Modelling the rms spectrum 

In the previous sections is was shown how the spectral variabil- 
ity is consistent with a model comprising a variable component 
V(-E) (the PLC; see Figure [Tsl and a constant component C{E) 
(the RDC; see Figure fTTl . Given these empirically derived spectral 
components it is possible to construct a model for the rms spectrum 
shown in Figure fT6l (see also Inoue & Matsumoto 2001). 

It is assumed that the spectrum is given by equation Q where 
both V{E) and C{E) include absorption effects. The (normalised) 
rms spectrum can be expressed as: 
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[F{E)] _ y^a[NyV(E)]^ + ^[NcC(E)]^ 



F{E) 



NvV{E) + NcC{E) 



(2) 



where a[F{E)] and F{E) represent the absolute rms ampli- 
tude and time-average of the spectrum over the observation, re- 
spectively. (Likewise for the two spectral components V{E) and 
C{E).) In terms of the model being discussed it is assumed that 
C{E) is not variable and so a[NcC{E)] = 0. Therefore: 



a[NvV{E)] 



NvV{E) + NcC{E) 



C{E)/F{E)) 



(3) 



and the term C {E) / F{E) is the fractional contribution of the con- 
stant component shown in Figure fTTI 

The above analysis assumes that the variable compo- 
nent changes only in normalisation and not in shape (i.e. 
a[NvV{E)]/NvV{E) = a[Nv]/N^), but is independent of the 
actual form of the spectrum V{E). The linearity of the flux-flux re- 
lations implies the variable component has a constant softness ratio 
with flux, and therefore its spectrum V{E) is independent of flux, 
validating this assumption. The first term on the right-hand side of 
equation|3]is thus simply a normalisation. Implicit in this analysis is 
the assumption that the absorption function is not time-variable. If 
the absorption function did vary this would introduce another term 
to the rms spectral model (Inoue & Matsumoto 2001), but such a 
term is not required by the data and thus it was assumed that the 
absorption did not vary. 

The histogram in Figure fT6l shows the resulting rms spectral 
model when the constant component derived in section l5!T1 (Fig- 
ure IIU is used. This very simple model clearly reproduces (to 
within ~ 1 per cent in JVar) the energy dependence of the variabil- 
ity amplitude on these timescales. The peak at ~ 1 keV is a result 
of the constant component C{E) being weakest at ~ 1 keV and the 
suppression at ~ 6.4 keV results from the strong iron line present 
in the constant component. The short timescale rms spectrum has a 
subtly different shape and thus it not so well reproduced using this 
simple model. However, it is known that at high frequencies the 
variability in different energy bands shows different PSD shapes 
and incoherent variability. Thus the short timescale rms spectrum 
will be affected by these additional effects, although their physi- 
cal origin is not clear. This could perhaps be explained if the RDC 
(or perhaps the warm absorber) varies on short timescales but these 
changes are 'washed out' over longer timescales. 
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Figure 17. The rm.s spectra, as shown in Fig. 1161 including only the low 
flux interval towards the end of the observation. 



longer timescale, shows a similar overall shape to the rms spec- 
trum from the entire observation (upper data from Figs. 1161 and 
I17t . Interestingly the prediction of the two component model dis- 
cussed above does not match these data so well. In particular the 
model predicts either too little rms around the iron line or too much 
in the soft X-ray band. This could be indicating that at low fluxes 
some variation in the shape/strength of the RDC were detected (see 
also Reynolds et al. 2003). However, this could also result from a 
random effect caused by the energy-dependent PSD (Vaughan et al. 
2003a,b). More observations at very low flux levels could confirm 
whether the RDC shows rapid variability. 



5.3.2 A low-flux rms spectrum 

The shape of the rms spectrum is clearly a function of timescale, as 
expected given the energy-dependent PSD (Vaughan et al. 2003a). 
It also shows subtle changes with time when the rms spectrum of 
each revolution of XMM-Newton data are compared (an rms spec- 
trum for the first revolution was shown in paper I). These may be 
caused by subtle changes in the RDC spectrum. As the relative 
strength of the RDC should be highest when the overall source flux 
is low, the low flux intervals of the light curve are the most promis- 
ing to search for changes in the RDC. 

The rms analysis described above was repeated using only 
the last 45 ks of data from rev 0303. During this time interval the 
source flux was rather low and in fact reached its minimum for 
the observation (Fig. I12t . The rms spectra from this interval are 
shown in Fig. 1171 The low-flux rms spectrum, particularly on the 



5.4 Time-resolved spectra 

For completeness the spectrum was examined in consecutive 10 ks 
time slices. Fig. llSl shows the spectra extracted from each of these 
time slices. Clearly the strength of the iron line (relative to the con- 
tinuum model) changes throughout the light curve. However, the 
only obvious, visible change is that described above; as the over- 
all X-ray flux decreases the iron line becomes more prominent be- 
cause it varies less than the continuum. For example, the spectrum 
for segment 303:1, which has one of the lowest fluxes, also has the 
most prominent line. 

It is possible that there are short-lived spectral features that 
would be 'washed out' in the time-averaged spectrum. Turner, 
Kraemer & Reeves (2003b) discuss the possible existence of tran- 
sient, redshifted iron lines. However, given that 32 independent 
spectra were examined here, the likelihood of an apparently signifi- 



16 Vaughan & Fabian 



cant feature appearing purely by chance is quite liigh. For example, 
the probability of detecting a feature, considered detected at 99.7 
per cent confidence in one individual spectrum, in one of the 32 
spectra is ~ 0.1 (using Pn f» 1 — (1 — Pi)^). By the same argu- 
ment only features found at > 3.5a significance in any individual 
spectrum should be considered as significant detections from the 
entire dataset (i.e. probability of false detection Pn < 0.01). There 
are no clear examples of strong, sharp features in these data except 
for the obvious iron line. The analyses of Fabian & Vaughan (2003) 
and Ballantyne et al. (2003) demonstrated that these time-resolved 
spectra can be fitted using relatively simple reflection models. In 
addition, there were no sharp features seen in the rms spectra (sec- 
tion |53j that would indicate a preferred energy for transient fea- 
tures. 



5.5 Principal Component Analysis 

The method of Principal Component Analysis (PCA) was used 
to try and isolate independently varying components present in 
the time-resolved spectra discussed above. PCA is a powerful sta- 
tistical tool used widely in the social sciences where it is often 
known as factor analysis. In practice PCA gives the eigenvalues 
and eigenvectors of the correlation matrix of the data. The data can 
be thought of as an m x n matrix comprising m observations of 
the n variables, in which case the correlation matrix would be a 
symmetric n x n matrix. 

The eigenvectors can be thought of as defining a new coor- 
dinate system, in the n-dimensional parameter space, which best 
describes the variance in the data. The first principal component, or 
PCI (the eigenvector with the highest eigenvalue), marks the direc- 
tion through the parameter space with the largest variance. The next 
Principal Component (PC2) marks the direction with the second 
largest amount of variance. The motivation behind PCA is to extract 
the (hopefully few) dominant correlations from a complex dataset. 
Francis & Wills (1999) provides a brief introduction to PCA as ap- 
plied to quasar spectra, while Whitney (1983) and Deeming (1964) 
illustrate more general applications of PCA in astronomy. 

If the variance in the spectra is caused by only a few indepen- 
dently varying spectral components then PCA should reveal them 
in its eigenvectors. The application of PCA to spectral variability 
of AGN was first discussed by Mittaz, Penston & Snijders (1990). 
The first few Principal Components (those representing most of the 
variance in the data) should reveal the shape of the relevant spectral 
components. The weaker Principal Components might be expected 
to be dominated by the photon noise in the spectra (which should 
be uncorrelated between each of the n energy bins in the m spectra 
and so should not distort the shape of the first Principal Compo- 
nents). 

This method was used on the m — 32 time-resolved spec- 
tra over the 3 — 10 keV range (binned such that they each spec- 
trum contained the same n = 21 energy bins). The spectra were 
unfolded using a F = 2 power-law model to convert them from 
counts to flux units (in EF{E) form). In the absence of sharp fea- 
tures in either the source spectrum or the detector response this 
should provide a reasonable estimate of the 'fluxed' data. The first 
three Principal Components are shown in Figure fT9l along with the 
mean and rms spectra derived from the same data. 

As can be clearly seen, PCI has a relatively flat spectrum and 
describes the vast majority of the variance in the data (96 per cent). 
The first Principal Component can therefore be identified with the 
variable PLC discussed above. All the remaining Principal Com- 
ponents each describe less than 1 per cent of the variance and most 
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Figure 19. Results of Principal Component Analysis (PCA). The top panel 
shows the mean spectrum (converted into EF{E) flux units). The next two 
panels show the rms and fractional rms of the data (after subtracting the 
contribution from photon noise). The bottom three panels show the spectra 
of the first three Principal Components and their contribution to the total 
variance of the data (Pvar). The solid line marks the zero level; points that 
lie on the same side of the zero line (all positive or all negative) are cor- 
related with one another while points that he on different sides of the zero 
line (differing signs) are anti-correlated. 



likely represent only the photon noise in the data. Thus the PCA 
confirms the above analyses and suggests the spectral variations 
in MCG~6-30-15 are a result of a single, varying continuum com- 
ponent (the PLC). Nearly identical results were found when the 
ranked data was used to produce the correlation matrix (which then 
represented the matrix of rank-order correlation coefficients). 



6 THE 2001-2000 DIFFERENCE SPECTRUM 

Figure 1201 shows the difference spectrum produced by subtract- 
ing the 3 — lOkeV EPIC pn:s spectrum taken in 2000 from the 
spectrum taken in 2001. A power-law provided a good fit (x^ — 
51.7/62 dof) with F = 2.18 ± 0.09 (compare with the results of 
section l5'.2.H . The only significant residuals appeared as an excess 
at ~ 6.6 keV. Including a narrow Gaussian improved the fit (x^ ~ 
44.1/60 dof), with a best-fitting energy E = 6.63 ± 0.12 keV. 
Including an absorption edge in the model did not significantly im- 
prove the upon this fit. Figure|^shows the ratio of the two spectra, 
confirming that the two spectra differ around E ~ 6.65 keV. 

An unresolved emission feature in the difference and ratio 
spectra could be due to a change in either the flux of the broad 
Fe emission line or the depth of the Fe Ka absorption. In order to 
produce an apparent emission feature in the first case only the blue 
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Figure 18. Top panel: EPIC pn 0.2 — 10 kcV light curve spanning all three XMM-Newton revolutions with 10 ks time intervals labelled. Bottom panels: EPIC 
pn spectra extracted from the 10 ks time intervals, shown as a ratio to a simple power-law model. The expected energy of the iron line (6.4 koV) is marked 
with a dotted line. 
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Figure 20. Difference spectrum produced by subtracting the 2000 data 
from the 2001 data. Only the EPIC pn:s spectrum is shown. The only sig- 
nificant residual (compared to the simple power-law continuum model) is 
an unresolved 'emission' feature at _E = 6.65 ± 0.05 keV. 
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Figure 21. Ratio of the EPIC pn:s spectral from the 2001 (upper data) and 
2000 (lower) observations. The top panel shows the raw spectra, the bottom 
panel shows their ratio. This again shows the unresolved 'emission' feature 
atE ^ 6.65 keV. 

wing of the emission line must be stronger in the 2001 observa- 
tion. In the second case the optical depth of the absorption must be 
deeper in the 2000 observation. An examination of the two spectra 
individually suggested the latter is more likely. A close examina- 
tion of changes in the RGS spectrum between the two observations 
would clarify this issue. 



7 REVIEW OF XMM-Newton RESULTS 

In this section the main results obtained from the long XMM- 
Newton observation of MCG-6-30-15 are summarised prior to the 
discussion in the following section. 

• The two XMM-Newton observations (taken in 2000 and 2001) 
sampled fairly typical 'states' of the source. The former observation 
sampled a period of lower flux than the latter. However, the long- 
term RXTE monitoring shows this can be attributed to short-term 
variability and does not imply a systematic difference between the 
states of the source during the two observations (section|3j- 

• The high energy spectrum obtained from BeppoSAX shows 
a strong Compton-reflection signature (paper I; Ballantyne et al. 
2003), as did the earlier XMM-NewtonlRXTE observation (WOl; 



Reynolds et al. 2003). There was no evidence for a low energy cut- 
off or roll-over in the continuum out to ~ 100 keV. 

• The RGS spectrum shows complex absorption by O VII, 
O VIII and Fe I as well as a range of other ions (Sako et al. 2003; 
Turner et al. 2003a). The opacity is concentrated mainly below 
~ 3 keV but still has an important effect on the 3 — 10 keV spec- 
trum (section l4.3.1> . See also Lee et al. (2001). 

• The fluorescent iron line is strong and broad (section l4!4l . The 
bulk of the line flux is resolved with EPIC. The emission peak con- 
centrated around 6.4 keV is resolved with a width FWHM ~ 
4.5 X 10* kms~^ , strongly indicating an origin within <J 100 rg. 
There is also a significant, asymmetric extension to lower energies 
that indicates strong gravitational redshifts. In addition there is a 
weak, intrinsically narrow core to the line emission (sec tions l4. 3.41 
andl4.4l see also Lee et al. 2002). 

• The best-fitting model for the 3 — 10 keV EPIC spectrum ex- 
plained the iron emission as from the surface of a relativistic accre- 
tion disc. The strong reflection explains the strength of both the iron 
line and the Compton reflection continuum. The best fitting model 
includes emission down to ~ 1-8 rg (section l4^ WOl; paper I). 

• The variations in X-ray luminosity show many striking simi- 
larities with those seen in GBHCs such as Cygnus X-1 (Vaughan et 
al. 2003a). In particular, the Power Spectral Density (PSD) function 
is similar to that expected by simply re-scaling the high/soft state 
PSD of Cyg X-1. The continuum variations are energy dependent 
and show the PSD is a function of energy and that the hard varia- 
tions are delayed with respect to the soft variations. Similar results 
have been found in other Seyfert Is (NGC 7469, Papadakis, Nandra 
& Kazanas 2001; NGC 4051, M'^Hardy et al. 2003). 

• Previous observations with A5C4 (Iwasawa et al. 1996, 1999; 
Shih et al. 2003) and RXTE (Lee et al. 2001; Vaughan & Edelson 
2001) showed the photon index of the 2 — lOkeV continuum to 
be correlated with its flux. The XMM-Newton observations confirm 
this and demonstrate the trend is reversed below ~ 1 keV, where 
the spectrum hardens with increasing flux (section l5!2l . The aver- 
age variability amplitude is highest in the range 1 — 2 keV, and 
is lowest at energies around the iron line (section l53l . 

• The variable spectrum can be decomposed into two compo- 
nents, a variable Power-law component (PLC; section l5'.2. It and a 
relatively constant Reflection Dominated Component (RDC; sec- 
tion lS.U . The spectral variability (at least on timescales ~ 10 ks) 
can be explained almost entirely by variations in the relative 
strength of these two components, caused solely by changes in the 
normalisation of the PLC (section [5.31 see also Shih et al. 2003 
and Fabian & Vaughan 2003). An analysis of the flux-flux plots 
(section ism see also Taylor et al. 2003) and the difference spectra 
(section l5.2.U shows the slope of the PLC remains approximately 
constant with flux. 

• The EPIC spectrum indicates there is resonance absorption by 
ionised Fe at ~ 6.7 keV (sections 14.3.21 and I4.4t . This was pre- 
dicted based on the presence of soft X-ray warm absorption (Matt 
1994; Sako et al. 2003) and has been observed in at least one 
other high quality EPIC spectrum (NGC 3783; Reeves et al. 2003). 
This resonance absorption appeared to vary between the two XMM- 
Newton observations (section|5l. 



8 DISCUSSION 

8.1 Fitting the iron line in MCG-6-30-15 

The broad iron line in MCG-6-30-15 has been unambiguously re- 
solved using XMM-NewtonfEPIC (WOl; paper I; this paper), con- 
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firming the previous results from ASCA (Tanaka et al. 1995; Iwa- 
sawa 1996), BeppoSAX (Guainazzi et al. 1999) and Chandra (Lee 
et al. 2002). However, deriving reliable line profile parameters is 
a considerable challenge even with the exceptionally high quality 
XMM-Newton+BeppoSAX data. From the whole of the 320 ks ob- 
servation the EPIC cameras collected about 6 x 10'' counts from 
the core of the line, peaking at 6.4 keV, and about three times as 
many from the broad red tail of the line emission. 

The foremost problem is that there are no regions of the X- 
ray spectrum unaffected by either the warm absorption (see sec- 
tion l4.3.l1 or the broad emission components. Thus it is not possi- 
ble to determine the underlying continuum without simultaneously 
constraining the other spectral components. Often used, simple 
techniques for determining the continuum by e.g. fitting a power- 
law over the 3 - 5 keV and 7 ~ 10 keV data (Nandra et al. 1997), 
will give a only crude first approximation of the continuum. In a 
source such as MCG-6-30-15, which possesses a strong warm ab- 
sorber, the absorption will cause the spectrum to curve even above 
~ 3 keV (see section U.B.ll and Fig.|3j. In the present paper this ef- 
fect was dealt with by including an absorption model derived from 
fits to the RGS data (section Pt.B.U and also allowing for additional 
absorption when fitting the EPIC data (section l431 . 

In addition to the distortion imposed on the continuum the 
other effect of warm absorption that has a significant impact on 
iron line studies is resonant line absorption at ~ 6.4 — 6.9 keV 
(section [4.3. 2> . The EPIC data indicate the presence of such ab- 
sorption; allowing for the Fe resonance absorption had a subtle but 
significant effect on the best-fitting emission line models (see sec- 
tion l4.4t . In particular, the line profile is slightly broader and bluer 
after accounting for the line absorption. 

Photoelectric absorption by Fe at energies ~ 7.1 — 9.3 keV 
could, in principle, also be caused by the warm absorbing gas 
and further confuse the continuum estimation (see e.g. Pounds & 
Reeves 2002). However, in the case of MCG-6-30-15 this is not a 
significant effect. The total opacity due to the Fe K edges, as pre- 
dicted based on fits to the soft X-ray RGS data, is negligible (sec- 
tion l4.3.2l . and spectral analysis of the EPIC data confirm this by 
ruling out the presence of deep Fe edges (section l4.5.H . The same 
is not true of all Seyfert Is galaxies. The columns of Fe ions derived 
using RGS data by Blustin et al. (2002) for NGC 3783, and by Sako 
et al. (2001) for IRAS 13349-1-2438 predict strong K edges with a 
total optical depth Tmax ~ 0.05 (see also Reeves et al. 2003). 

Further complications include the presence of the reflection 
continuum (section l4.3.3> and possible weak, narrow components 
to the line emission (section I4.3.4t . The approach used in the 
present paper was an attempt not to 'unambiguously isolate' the 
iron line but to model it after allowing for these complicating fac- 
tors. In order to limit the range of models and free parameters, 
constraints obtained from other instruments (such as the RGS and 
BeppoSAX spectra) were used where possible. The result of this 
spectral fitting confirms the previous analyses in WOl and paper I: 
The best-fitting emission model comprises a power-law continuum 
plus strong reflection from a weakly ionised, relativistic accretion 
disc fsection l431 . This model requires the disc emission to be cen- 
trally concentrated and extend in to ~ 2rg. Attempts to 'explain 
away' the red wing of the line using enhanced absorption or bro- 
ken continua (section |431 did not substantially alter the main re- 
sults. Additionally, alternative emission components such as iron 
line blends did not fit the data (section l44l while emission from 
soft X-ray lines or blackbodies led to physically implausible results 
(section l4.7t . 

The only non-standard model that could fit the data with plau- 



sible parameters was the partial covering model (section l4?7l . How- 
ever, as shown also by Reynolds et al. (2003), this model is strongly 
at odds with the higher energy data. Thus, a partially covered con- 
tinuum does not provide a satisfactory explanation of the X-ray 
spectrum of MCG-6-30-15. However, it is possible that the source 
does contains some partially covered emission in addition to the 
disc reflection. Fitting the data with a relativistic disc and also al- 
lowing for partial covering suggested the possible impact of the 
partial covering is minimal (section l477l . 

8.2 Spectral variability of MCG-6-30-15 

The spectral variability of the source can be explained in terms of a 
two-component model (section |5!T1 see also M'^Hardy et al. 1998; 
Shih et al. 2002; Taylor et al. 2003; Fabian & Vaughan 2003). This 
is the simplest model consistent with the flux-flux analysis. 

In this model the two spectral components are a power-law 
component (PLC) and a reflection dominated component (RDC) 
which carries most of the iron line flux. Both of these emission 
spectra are viewed through the warm absorbing gas. Flux variations 
are caused primarily by changes in the normalisation of the PLC, 
with the RDC remaining relatively constant. Such a model explains 
many aspects of the spectral variability including the linearity of 
the flux-flux plots (section [5. 1> . the rms spectrum (section f5.3t . 
the correlation between 2 — 10 keV spectral slope and flux (sec- 
tion l5.2> and the lack of strong iron line variations (Reynolds 2000; 
Vaughan & Edelson 2001; paper I). The relative constancy of the 
RDC causes variations above ~ 2 keV (and particularly around the 
iron line) to be suppressed, and also variations below ~ 1 keV are 
increasingly suppressed (section [5.31 . Fabian & Vaughan (2003) 
previously applied this model to direct fitting of the time-resolved 
EPIC spectra, where it gave a quite reasonable description of the 
time varying spectrum. 

The spectrum of the PLC was uncovered using the hi-low dif- 
ference spectrum (Fabian & Vaughan; section \5l2\ and the spec- 
trum of the RDC was obtained by isolating the constant offset of 
the flux-flux plots (Taylor et al. 2003; section lsTll . This latter tech- 
nique revealed the spectrum of the constant RDC independent of 
any spectral fitting, simply by analysing the light curves using flux- 
flux plots. The result clearly revealed a spectral form strongly sug- 
gestive of reflection, with a prominent iron line (as suggested by 
Fabian & Vaughan 2003). The RDC also shows a 'soft excess' 
below ~ 1 keV which could also be emission from the reflect- 
ing disc if it is weakly ionised. The model derived from fitting the 
3—10 keV EPIC spectra is in rough qualitative agreement with 
the RDC spectrum obtained independently (see Fie. l22t . Thus the 
model derived from spectral fitting (section |4} is in broad agree- 
ment with the model invoked to explain the spectral variability (sec- 
tion|5}. 



8.2.1 Variability constraints on complex absorption models 

The warm absorption spectrum does not appear to vary substan- 
tially (specifically, the optical depths of the absorption features re- 
main approximately constant). For example, ratio plots of 'high' 
and 'low' spectral (section ls!2l show smooth profiles over the deep- 
est absorption features, as does the rms spectrum (section|53}- This 
result is consistent with the studies of the deep warm absorber in 
NGC 3783 using Chandra (Netzer et al. 2003) and XMM-Newton 
(Behar et al. 2003) which showed the X-ray absorption remained 
unchanged despite large changes in source luminosity. 
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Figure 22. Fractional contribution to tlie spectrum of tlie constant compo- 
nent deduced as in 

Fig.ini The dotted curve shows the prediction based on 
the model deiived from fitting the EPIC spectra in the 3 — 10 keV range 
only. Clearly the model has qualitatively the same shape, although the de- 
tails do not match, particularly at lower energies. However, given that the 
model was derived without using the spectral data below 3 kcV, and there 
is an intrinsic uncertainty in the normalisation of the RDC (section im . the 
agreement is rather interesting. 

Furthermore, the spectral variability constrained the possible 
effect of warm or partial absorption on the iron line. Neither of 
the two spectral components isolated through spectral variability 
analysis resemble the absorbed part of the continuum in the par- 
tial covering model. The depth of any line-of-sight Fe K edge was 
constrained by the difference spectrum analysis (section l5!2l to be 
r < 0.02 at 7.1 keV tor < 0.07 at 8.0 keV. Fitting the partial cov- 
ering model to the various difference spectra gave constraints com- 
parable to those based on the time-averaged spectrum (section l477l . 
As discussed above the effect of including such additional absorp- 
tion (either warm or partial) on the derived iron line parameters was 
negligible. 

Alternative models for the spectral variability of MCG-6- 
30-15 in which the variations are caused primarily by changes 
in the warm absorber (Nandra et al. 1990; Inoue & Matsumoto 
2003) are in stark contrast with the data. This model predicts a 
power-law flux-flux relation (since the flux is given by f{E, t) = 
S{E) exp(— r(i?, t)), where S{E) is the underlying the emission 
spectrum and t{E, t) is the time-variable absorption optical depth). 
Such a model would also predict strong features in the rms spec- 
trum corresponding to the deepest edges in the absorption spectrum 
(paper 1). Both of these were ruled out. 

8.3 Outstanding issues with the disc model 

8.3.1 Lack of iron line variations 

Although the model outlined above explains many aspects of the 
EPIC spectrum and its variability, there are some serious out- 
standing problems. The most significant of these is how the RDC 
flux, including the iron line, remains so steady in the presence 
of large changes in the flux of the PLC. In standard accretion 
disc/reflection models the reflection spectrum (including iron line) 
is fundamentally driven by the PLC luminosity. Thus, since the 
RDC is thought to originate close to the SMBH (and the source 
of the PLC) the fluxes of the two components should be tightly 
coupled. The lack of correlations between the iron line and con- 



tinuum is a long-standing issue (Iwasawa et al. 1996, 1999; Lee et 
al. 2000; Reynolds 2000; Vaughan & Edelson 2001). These XMM- 
Newton observations demonstrate more clearly than before that not 
just the line but the bulk of the reflected emission (the RDC) shows 
little short-timescale variability. The physics responsible for the re- 
duced variability of the RDC is unclear. 

During the 2001 XMM-Newton the total EPIC count rate from 
the iron line was '--^ 0.1 ct s^^, with ^ 0.02 ct s^^ coming from 
the resolved, quasi-Gaussian core on the line and the rest coming 
from the broad, red wing. As a result of this low count rate, com- 
pounded by the modelling difficulties, it is difficult to place firm 
limits on weak and/or extremely rapid (faster than ~ 10'* s) line 
variations. Intriguingly, Reynolds et al. (2003) claim that variations 
in the line/RDC become apparent only when the source flux is very 
low (see also Iwasawa et al. 1996). 

One possibility is the progressive ionisation of disc surface, as 
discussed by e.g. Reynolds (2000), Nayakshin & Kazanas (2002) 
and Ballantyne & Ross (2002). In this scenario an increase in the 
ionising luminosity leads to little change in the flux of the line be- 
cause more of the Fe ions in the surface layers of the disc become 
fully-stripped of electrons. Under some circumstances this model 
can explain the lack of iron line variations (Reynolds 2000). How- 
ever, the spectral fitting analysis (section|4} suggests the emitting 
region of the disc is only weakly ionised. Time-resolved spectral fit- 
ting of the XMM-Newton data using an ionised disc model (Ballan- 
tyne et al. 2003) showed that ionisation effects could not reproduce 
the constancy of the RDC. A second reflector that was intrinsically 
constant was also required. Thus, ionisation effects alone are insuf- 
ficient to explain the lack of iron line variations. 

One intriguing alternative possibility is that gravitational light 
bending, which may be strong if the disc extends to ~ 2 rg, is 
(partially) responsible (Fabian & Vaughan 2003). See Miniutti & 
Fabian (2003) for a detailed discussion of this model. Such a model 
may naturally explain the enhanced reflection strength (R ^ 2; sec- 
tion l4.3.3l and iron line. 

An interesting and related question is: how does the PLC re- 
tain its spectral shape whilst undergoing factors of 5 changes in 
luminosity? Both direct time-resolved spectral fitting (e.g. Fabian 
& Vaughan 2003) and flux-flux analyses ( section IsTTl Taylor et al. 
2003) show the slope of the PLC to stay remarkably constant. Simi- 
lar results have been seen in ultrasoft Seyfert galaxies (e.g. Vaughan 
et al. 2002) where the primary power-law slope stays constant de- 
spite large amplitude luminosity variations. This poses an interest- 
ing challenge to models of the origin of the primary X-rays which 
typically predict some luminosity-correlated variation in slope (e.g. 
Zdziarski et al. 2003). The observations are suggestive of some 
kind of 'Compton thermostat' (Haardt, Marachi & Ghisellini 1994; 
Pietrini & Krolik 1995). Alternatively, if the apparent variations 
in the PLC flux are due to geometrical effects near the SMBH, 
the PLC might be expected to vary achromatically (Miniutti et al. 
2003). 

8.3.2 Inner edge of the disc 

The intricacies of modelling the complex spectrum of MCG-6-30- 
15 have made it difficult to tap the potential of the iron line pro- 
file for probing the region of strong gravity around the SMBH. 
An often asked question is: how far in towards the SMBH does 
the disc extend? The answer to such a question may reveal the 
spin of the SMBH (e.g. Stella 1990; Martocchia et al. 2002). The 
strength/shape of the red wing of the line is determined by the emis- 
sivity of the innermost disc and so are the primary diagnostics of 
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the position of the disc's inner boundary. The problem is that the 
very broad, low-contrast tail to the line is difficult to discern from 
the spectral curvature caused by the warm absorber (see Fig.|3j- 
Thus it is not possible to unambiguously determine the inner disc 
edge (unless the transmission of the warm absorber is known to 
very high accuracy). For example, the method of Bromley, Miller 
& Pariv (1998) works by identifying the 'minimum energy' of the 
red wing of the redshifted line emission. However, no such energy 
can be identified without detailed, simultaneous modelling of the 
absorption and the reflection continuum. 

This does not mean that all is lost. The method of fitting the re- 
flection emission while simultaneously fitting/constraining the ab- 
sorbed continuum did require emission in to ~ 2 rg (WOl; paper I; 
section l431 . Importantly, this result was robust to various possible 
effects that might bias the continuum determination (sections 14.51 
and l4.7t . This is strong but not conclusive evidence that the SMBH 
is spinning. An alternative approach to fitting the red wing of the 
line (perhaps isolating it through variability if the red wing shows 
variations on very short timescales not yet probed) would bolster 
the determination of its profile. 

The second and more fundamental objection is that some 
emission is possible from within 6rg even if the SMBH is not 
spinning (see Reynolds & Begelman 1997; Young, Ross & Fabian 
1998; Krolik & Hawley 2002; Merloni & Fabian 2003). Further 
refinement of the theory combined with a further increase in obser- 
vational sensitivity should be able to settle this issue (see Young & 
Reynolds 2000; Yaqoob 2001). 



of broad lines in many Seyfert Is is therefore quite puzzling and 
requires that either there simply is not enough cool matter extend- 
ing close to the SMBH to produce the line or that the line photons 
are somehow lost (absorbed/scattered) on their way out of the in- 
ner regions. See also Reeves et al. (2003) for a discussion of this 
problem. 

It is important to note that there could be some observational 
biases present. By their very nature highly broadened lines are dif- 
ficult to isolate (as discussed above). The broad line in MCG-6- 
30-15 is perhaps the most clear example in part due to the strength 
of the reflection spectrum. The high level of reflection (R > 2) 
means that the equivalent width of the broad line is rather high 
{EW ~ 400 eV). Other Seyfert galaxies with more typical re- 
flection strengths (R ~ 1) would be expected to have comparably 
weaker broad lines {EW ~ 150 eV; George & Fabian 1991). Even 
with the high throughput of the EPIC cameras on-board XMM- 
Newton, a typical length observation (~ 40 ks) of a bright Seyfert 1 
might not clearly reveal a redshifted line with EW ^ 150 eV. This 
is especially true if the continuum is absorbed. Thus, most observa- 
tions of Seyfert Is are unlikely to show broad lines; the lack of an 
obvious broad line in the spectrum does not rule out its existence. 
These observations can constrain the strength of highly broadened 
lines (e.g. upper limits on the equivalent width of a Laor profile) 
and so place useful limits on the line emission. This will help ad- 
dress the question of whether broad lines are a rare or common 
occurrence. 



8.3.3 Why do only some Seyferts show broad lines? 

Recent, high-quality Chandra and XMM-Newton observations of 
bright Seyfert galaxies seem to have produced a mixed bag of 
iron lines. Examples of strong, narrow lines have been found in 
NGC 3783 (Kaspi et al. 2002; Reeves et al. 2003), NGC 5548 
(Yaqoob et al. 2001; Pounds et al. 2003a), NGC 4151 (Schurch 
et al. 2003), and NGC 3227 (Reeves 2003). Examples of highly 
broadened iron lines have been seen in MCG-5-23-16 (Dewangan 
et al. 2003), NGC 3516 (Turner et al. 2003c), IRAS 18325-5926 
(Iwasawa et al. in prep.), Q 0056-363 (Porquet & Reeves 2003), 
IRAS 13349-1-2438 (Longinotti et al. 2003) and Mrk 766 (Mason 
et al. 2003; Pounds et al. 2003b). The line in MCG-6-30-15 is 
the broadest of the known broad lines. There are some objects for 
which the line profile is in dispute (e.g. Mrk 205 and Mrk 509, 
see Pounds et al. 2001; Reeves et al. 2001; Page et al. 2003c). 
Even more confusing is the possible detection of rapidly variable 
yet narrow iron lines in Mrk 841 (Petrucci et al. 2002) and NGC 
7314 (Yaqoob et al. 2003). Recent studies of GBHCs have also re- 
vealed relativistically broadened iron lines in a few cases (Miller et 
al. 2003; Miller et al. 2002). 

This raises two obvious, related questions. The first is: why 
do some Seyferts show strong, broad lines while others show only 
narrow lines (with any highly broadened component being weak or 
absent)? The second question is: how can some Seyferts not have 
broad lines? The high luminosities of AGN in general argues that 
there must be a large amount of mass flowing deep into the po- 
tential well of the SMBH. Furthermore, the rapid X-ray variability 
argues that the X-ray source is spatially compact and probably also 
located close to the SMBH. There should be a substantial amount 
of relatively cool accreting matter in the region of strong gravity 
around the SMBH, close to the primary X-ray source. These are 
the only two ingredients thought necessary to produce a broad iron 
line (Fabian et al. 1989; Reynolds & Nowak 2003). The absence 



8.4 Iron resonance absorption 

The fits to the iron line allowing for resonant absorption suggest 
the presence of an additional absorber whose dominant effect is 
Fe XXV Hea absorption at ~ 6.7 keV. This absorption appears 
to have varied between the 2000 and 2001 XMM-Newton obser- 
vations (section |6}. Such an absorption system can be modelled 
(using the CLOUDY grids of Turner et al. 2003a) with a high 
column (TVh > 10^^ cm~^ ) and high ionisation (log(^) « 3) 
warm absorber. This is similar to the results obtained for another 
bright Seyfert 1, NGC 3783 (Reeves et al. 2003) except that this 
very high ionisation absorption does not significantly affect the 
low energy spectrum. Similar features may have been observed in 
other Seyferts (NGC 3516, Nandra et al. 1999; IRAS 13349-2438, 
Longinotti et al. 2003). 

Another approach is to estimate the column density using just 
the equivalent width (EW ~ 10 eV) of the absorption line and 
the assumption that He-like Fe is the dominant ion in the absorbing 
material. Using the oscillator strength of the Fe XXV Hea line from 
Nahar & Pradhan (1999) the column density can be estimated as- 
suming the line is unsaturated (linear part of the 'curve of growth' ; 
Spitzer 1978). This gives Np^-i-24, ~ lO^'^ cm~^ . If the line is satu- 
rated the column density will be higher The predicted depth of the 
corresponding Fe XXV absorption edge (at 8.83 keV) is r ~ 0.05, 
consistent with the limits obtained from the data (section l4.5.1> . As- 
suming solar abundances and that half all the Fe ions are in the He- 
like state this column corresponds to A^'h > 4 x lO^'^ cm~^ . This 
compares reasonably to the estimate derived from the CLOUDY 
model, but should really be considered only a lower limit on the 
total column density since it does not account for other (lower os- 
cillator strength) lines, possible line saturation or emission line fill- 
ing. 
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8.5 The Future 

The long XMM-Newton observation of MCG-6-30-15 has yielded 
a vast amount of detailed information (on absorption and emission 
spectra, variability timescales and delays and spectral changes). 
Only by making full use of the long and simultaneous observa- 
tions (CCD and grating X-ray spectra) possible with XMM-Newton 
(and simultaneous BeppoSAX observations) could all of these be 
achieved. Comparable observations have been performed/planned 
for a small sample of other Seyferts (for example: NGC 3783, De- 
bar et al. 2003; Reeves et al. 2003) which will shed light on the 
detailed X-ray properties of Seyferts in general. Complimentary to 
these are the Chandra long-looks that have been particularly en- 
lightening for warm absorber studies (e.g. Kaspi et al. 2002). Such 
observations also demonstrate the limits of present instrumentation. 
For example, without higher resolution, high S/N spectroscopy it is 
difficult to realise the diagnostic potential of the 6.4 — 6.9 keV iron 
resonance absorption lines. 

The presence of such a large column of highly ionised iron 
makes MCG-6-30-15 a good choice for high resolution Fe-band 
spectroscopy with the XRS on-board Astro-E 2. This has the capa- 
bility (with resolution Ai5 ~ 7 eV and good throughput) to mea- 
sure the depths and velocity shifts of some of the individual res- 
onance lines. Indeed, the combination of the high resolution XRS 
spectrum of the iron absorption/emission, the high energy HXD 
spectrum which will constrain the continuum and reflection, and 
the high throughput XIS spectrum should provide the next major 
advance towards making full use of the information contained in 
the X-ray spectra of Seyferts. 

Aside from future missions, the present data indicate the most 
interesting observations are often made when the source is very 
faint. The spectral analysis presented above showed that when 
MCG-6-30-15 is faintest, the iron line is strongest relative to the 
continuum, making it a much higher contrast feature. The rare ob- 
servations that show apparent changes in iron line properties tend 
to occur during periods of unusually low flux (Iwasawa et al. 1996; 
Reynolds et al. 2003). Observations of NGC 4051 during its pro- 
longed 'low states' have revealed other interesting phenomena, 
such as the warm absorbing gas in emission (Uttley et al. 2003). 
Thus it would seem that XMM-Newton and Chandra observations 
of bright Seyfert Is, taken during periods of unusually low activity, 
may prove particularly revealing. 



9 CONCLUSIONS 

The long XMM-Newton observation of MCG-6-30-15 has con- 
firmed the existence of a highly broadened iron line. The main re- 
sults of the present investigation are as follows: 

(i) It is not yet possible, given the present state of knowledge 
of the detailed properties of the warm absorber, to unambiguously 
dis-entangle the red wing of the line from the absorbed continuum. 

(ii) Nevertheless, the best-fitting model requires the reflec- 
tion/line emitting region to extend inwards to ~ 2rg. This result 
was robust to various possible biases in the continuum determina- 
tion. 

(iii) Alternative models for the spectral region around the iron 
line, not allowing for the effects of strong gravity, gave unsatisfac- 
tory results. 

(iv) The spectral variability of MCG-6-30-15, including the cor- 
relation between 2 — lOkeV slope and flux and the reduced vari- 



ability of the iron line, has been explained using a two-component 
model. 

(v) The relative constancy of the reflection dominated compo- 
nent (RDC), in the presence of large variations in the power-law 
component (PLC), is the primary cause of the spectral variations. 
The reasons for the lack of RDC variations are unclear. A model in- 
voking gravitational light bending near the SMBH can qualitatively 
explain the suppressed variations, the relative strength of the RDC 
and the small inner disc radius. 
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Figure Al. Confidence contours from fitting the 3C 273 EPIC spectra 
(from revs. 0094 — 0096) over the 3 — 10 keV range with a power-law 
model. 




Channel energy (keV) 

Figure A2. 3-10 keV EPIC spectra of 3C 273 from revs. 094 - 096. The 
four sets of data are for pn single (filled circles) and double events (filled 
triangles), MOSl (hollow circles) and MOS2 (hollow triangles) The spectra 
were fitted with a model comprising a power-law and a Gaussian emission 
line. 



APPENDIX A: ANALYSIS OF EPIC CALIBRATION DATA 

This section describes an analysis of EPIC data taken from obser- 
vations of 3C 273 (z — 0.158). This source is bright, point-like 
(compared to the EPIC PSF) and shows a relatively simple energy 
spectrum in the EPIC bandpass (i.e. approximately power-law with 
no sharp spectral features), making it useful for assessing the accu- 
racy of the EPIC spectral calibration. 

Five observations of 3C 273 were analysed. These observa- 
tions were taken primarily for performance verification and cali- 
bration purposes and the observation details are given in TablelAll 
During all five observations the source was positioned on-axis and 
all three EPIC cameras used small window mode and medium fil- 
ter (i.e. the same instrumental set-up as the MCG-6-30-15 observa- 
tions). The only exception is the first half of the rev. 277 observa- 
tion during which the MOS cameras used the thin filter (data from 
this part of the observation were ignored). This later observation 
was used by Molendi & Sembay (2003) to assess the EPIC calibra- 
tion. 

Source and background spectra were extracted exactly as for 
the MCG-6-30-15 observation. Only single pixel event were used 
from the MOS to reduce the effect of pile-up in this bright source 



(see Molendi & Sembay 2003). The grouped spectra were then fit- 
ted in the 3 — 10 keV range with a simple power-law modified by 
Galactic absorption (A^h = 1.79 x 10^° cm"^ ; Dickey & Lock- 
man 1990). The power-law indices and normalisations were left 
independent for each camera. In all five observations the spectra 
from each camera showed weak but positive residuals at 5 — 6 keV, 
suggesting the presence of iron line emission. Previous studies of 
3C 273 (e.g. Turner et al. 1990; Grandi et al. 1997; Reeves 2003) 
also suggested the presence of weak iron emission. A Gaussian was 
therefore included in the model at 6.4 keV (in the source frame) to 
account for iron emission (with equivalent width EW ~ 20 eV). 
This simple model provided an acceptable fit to all five sets of spec- 
tra (i.e. xi < 1-0) and showed that the four sets of spectra taken 
during rev. 094 — 096 were very similar (i.e. the source showed 
no spectral variability during this period). Therefore the data from 
the four observations spanning these revolutions were combined to 
produce a high signal-to-noise spectrum for each EPIC camera. 

Fig. lAll shows a comparison of the photon indices and normal- 
isations for the three EPIC spectra from the rev. 094 — 096 data. 
Similar results were obtained from the rev. 277 data. Clearly the 
spectral indices are discrepant, with the MOS 1 spectral index sig- 
nificantly lower than M0S2, which in turn is lower than the pn. The 
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Table Al. Details of observations used for calibration purposes. 



Observation Observing time Duration 
Source Revolution ID (UT) (s) 



3C 273 


094 


0126700301 


2000 Jun 13 23:39:53- 


-2000 Jun 14 20:05:49 


73556 


3C 273 


095 


0126700601 


2000 Jun 15 12:58:18- 


-21:35:30 


31032 


3C 273 


095 


0126700701 


2000 Jun 15 23:32:02- 


-2000 Jun 16 09:37:48 


36346 


3C 273 


096 


0126700801 


2000 Jun 17 23:24:14- 


-2000 Jun 18 19:50:15 


73561 


3C 273 


277 


0136550101 


2001 Jun 13 07:14:26- 


-2001 Jun 14 08:10:31 


89765 



flatter index derived from MOSl is a known calibration problem 
(Molendi & Sembay 2003), and the other three spectra are within 
Ar « 0.07 of each other. As there is likely to be some spectral 
distortion due to pile-up in the MOS data, even in the single pixel 
events, the differences in slope between the MOS cameras and the 
pn might be expected to be slightly smaller for a fainter (i.e. less 
piled-up) source such as MCG-6-30-15. Fie. IA2I shows the data 
and residuals from the power-law plus Gaussian fit. Clearly this 
simple model reproduced the EPIC spectra to better than ±5 per 
cent over the 3 — 10 keV range. Excepting these systematic differ- 
ences in spectral slope between the three cameras, the EPIC spec- 
tral calibration appears good to ~ 5 per cent in the 3 — 10 keV band 
(Fig. lA2t . The only noticeable sharp residual in the 3 — 5.5 keV 
range occurs at « 4 keV in the observed frame. A shallow absorp- 
tion edge fits this feature (E = 4.UtaT7 keV, r = 0.023^0:0^5). 
This is possibly an instrumental feature caused by Ca deposits on 
the mirror surfaces, but is only present as a ~ 3 per cent effect. 

Extrapolating the above model down to lower energies re- 
veals and upturn in the data, i.e. 3C 273 has a 'soft excess'. The 
0.5 — 3kcV spectrum can be fitted with a broken power-law or 
two blackbody components to account for the curved continuum. 
The parameters were allowed to vary independently for each cam- 
era and again the spectral slopes were significantly different. The 
fit was poor but the residuals are generally confined to the ±4 per 
cent level. The most noticeable exceptions are near the instrumen- 
tal edges (O K-edge at 0.538 keV, Si K-edge at 1.846 keV and Au 
M-edge at 2.21 keV), particularly evident in the high S/N pn spec- 
tra. 

This paper has been typeset from a TgX/ ET^ file prepared by the 
author. 



